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Abstract 

We investigate the long time behavior of the critical mass Patlak-Keller-Segel equation. This equation has a 
one parameter family of steady-state solutions q\, X > 0, with thick tails whose second moment is not bounded. 
We show that these steady state solutions are stable, and find basins of attraction for them using an entropy 
functional Hx coming from the critical fast diffusion equation in R 2 . We construct solutions of Patlak-Keller- 
Segel equation satisfying an entropy-entropy dissipation inequality for H\. While the entropy dissipation for 
H\ is strictly positive, it turns out to be a difference of two terms, neither of which need to be small when the 
dissipation is small. We introduce a strategy of controlled concentration to deal with this issue, and then use the 
regularity obtained from the entropy-entropy dissipation inequality to prove the existence of basins of attraction 
for each stationary state composed by certain initial data converging towards q\. In the present paper, we do 
not provide any estimate of the rate of convergence, but we discuss how this would result from a stability result 
for a certain sharp Gagliardo-Nirenberg-Sobolev inequality. 

Mathematics subject classification numbers: 15A45, 49M20 

1 Introduction 

1.1 The PKS system and its critical mass 

The Patlak-Keller-Segel system HZ] is one of the simplest models of chemotaxis, describing the evolution of the 
population density of a cell colony which is diffusing across a two dimensional surface. In addition to the diffusion, 
as the cells move across the surface, they continually emit a chemical attractant, which itself diffuses across the 
surface. The cells tend to move towards higher concentrations of the attractant, and this induces a drift term 
tending to concentrate the population, and countering the spreading effects of the diffusion. A model organism for 
this type of behavior is the Dictyostelium Discoideum which segregates cyclic adenosine monophosphate, another 
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important example of chemotactic movement are endothelial cells who react to VEGF to form blood vessels. See 
[24l |40] for recent reviews on chemotaxis models. 

The model is mathematically interesting on account of this competition between the concentrating effects of the 
drift induced by the chemical attractant and the spreading effects of the diffusion, and there is a critical value to the 
total mass of the initial data, so that for masses above this value, the concentration wins, and the density collapses 
in a finite time. However, for masses below this critical mass, diffusion dominates, and the colony smoothly diffuses 
off to infinity. At the critical mass, there is a continuous family of stationary solutions, and this paper is concerned 
with determining their stability properties, and since they all turn out to be stable, basins of attraction for each 
of them. We begin by introducing the model and the critical mass associated with it. If p denotes the population 
density, and c the concentration of the chemical attractant, the system of equations is 



^(t,x) = div [Vp(t,x 



c(t, x) - 



1 

2^ 



/II 2 
p (x) > 



p{t, x)Wc(t, x)] 
log \x - y\p(t,y) dy , 



t > 0, x € 
t > , x e 
x e R 2 . 



(1.1) 



with an appropriate choices of units, so that all dimensional constants are unity. 

In most of this paper, we consider initial data po that belongs to L 1 (R 2 , log(e+ \x\ 2 ) dx), and such that p logpo 
is integrable. The relevance of these conditions shall be explained shortly, but at the very least, they insure that 
c(0, x) is well defined. It will sometimes be convenient to write the second equation in in the compact form 
c(t, x) — G * p(t, x) where G(x) = — 1/(2 tt) log \x\ is the Green's function for —A in R 2 . That is, —Ac = p. 

Also throughout the paper, the term density shall always refer to a non- negative integrable function on R 2 , and 
we shall use the term mass to refer to the total integral of a density p. Because of the divergence form structure of 
the system, solutions formally satisfy the conservation of mass 



p(t, x) dx 



p {x) dx := M 



for all t > 0; i.e., the mass M is conserved in time. 

The PKS system can be rewritten advantageously as follows: Introduce the free energy functional J-pks 

•FpKsH = / Pi x ) p{ x ) dx + -— / / p(x)\og\x-y\p(y)dxdy . 

The first integral is well defined if plogp is integrable, and the positive part of p(x)\og\x — y\p{y) is integrable 
when p belongs to i 1 (R 2 , log(e + \x\ 2 ) dx), so that the second integral is at least well-defined under this condition. 

Now suppose that the density p belongs to L 1 (R 2 ,log(e + \x\ 2 ) dx), and moreover, plogp is integrable. Then a 
simple formal calculation shows that for all u £ C^°(R 2 ) with zero mean, 



lim - (.FpksIp + eu \ - -^pksIp]) = / 



dp 



(x) u{x) dx 



where 



fr^PKS P ^ . = \ ogp {x) + / lot; ./■ - u i>{y)dy = l>;xp(.r) - ,>(.r) . 



It is then easy to see that the evolution equation in (11.1[) can be rewritten as 



Sp 



It follows that at least along well-behaved classical solutions (for which we may integrate by parts), 

2 



dt 



•^PKS [/<>(*)] 



p(t, x) 



_^PKSPWL s 

V (x) 

dp 



dx . 



(1.2) 



(1.3) 
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In particular, along such solutions, t H> J-pKs[p{t)] is monotone non-increasing. The key to exploiting this mono- 
tonicity, as discovered in [21], is the sharp logarithmic Hardy-Littlewood-Sobolev (Log HLS) inequality [31 [12]: 

1.1 LEMMA (Logarithmic Hardy-Littlewood-Sobolev inequality). Let f be a non-negative measurable function 
in R 2 such that /log/ and /log(e + \x\ 2 ) belong to L 1 (M 2 ). Let M := / R2 / dx. Then 

f(x) log f(x) dx + A // /(a) log |x - y|/(2/) dx dy > - C(M) , (1.4) 

with C{M) := M (1 +log7r — log(ilf)). There is equality if and only if /(x) = £>a(£ — £o) /or some A > and some 
xo € R 2 , iw/iere 

M A 

0\{X) := — — — —2 • 

7T (A+ |x| 2 ) 2 

Following [5T], one may apply sharp log HLS inequality (|1.4p to deduce that 



^KB[p] = g(/ R3 P(.) I0gp(x)dx + A|| R2 



p(x) log \x - y\p(y) dx dy 



>-M f 1 _M\ I p{x)logp{x)dx _ (15) 

8tt V 8nJ J R 2 

It follows from this and the monotonicity of J-pKs[p(t)] that for solutions p of the PKS system for which M < 8ir, 

^r/M f / \ ^ { \ i 87rJ"[p ] -MC(M) 
f p * := / P(t,x )\og p(t,x dx< ^ v ; . 

Therefore, for M < 8ir, the entropy £ [p(t)] stays bounded from above, uniformly in time. This precludes the 
collapse of mass into a point mass for such initial data. In |21i [7], this formal analysis is made rigorous, and the 
global existence of solutions below the critical value 8-7T and a number of their properties as well are established. 

Previous work in this direction, by Jager and Luckhaus [25) , had shown that for initial data of sufficiently small 
mass, the entropy £[p(t)] stayed bounded uniformly in t. Their analysis used the Gagliardo-Nirenberg-Sobolev 
inequality for functions / in R 2 that bounds ||/||4 in terms of ||V/||2 and ||/||2, and not the Log HLS inequality, 
but their global existence result requires the mass to lie below a threshold that is strictly less than 8-7T. 

That 8n is the actual critical value at which diffusive and concentrating effects arc balanced, and not only a 
better lower bound, can be seen by computing moments: When the initial data has a finite second moment, and 
M > 8tt such collapse, or "blow-up" docs indeed occur in a finite time. To see this, we first note a weak formulation 
of our the PKS evolution equation that will useful to us later on. Let tp be any test function. Then 



d_ 

dt 



ip(x)p(t, x) dx — / Aip(x)p(t, x) dx 



1 f , (Vy>(a) - Vy>(|/ )) ■ (a - y) 

47t/r2 x r2 \x-y 



P(t,x) — — 2 p(t,y)dxdy. (1.6) 



In addition to the usual integration by parts, we have symmetrized the second term on the right in x and y. 

d f 

Fixing any a 6 R 2 and taking ipix) = a ■ x, we see from (ll.6[) that — / xp(t,x) dx = 0; i.e., the center of mass 

dt J R 2 

is conserved. Due to the translational invariance, we henceforth assume zero center of mass. More interestingly, 
taking ip(x) = |x| 2 , so that (Vip(x) — Vip(y)) ■ (x — y) = 2, we find 



1 / \x\ 2 p(t, x) dx = AM ~ -Uf 2 = 4M ( 1 - ^ 
dt J R 2 2vr V 8?r 
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Thus, if M > 8n, the right hand side is strictly negative, and this shows that the second moment of p(t) reaches 
zero in a finite time if initially bounded, or else some sort of singularity develops that would invalidate the formal 
calculation we have just made. 

Thus, the mass value M = 8ir is the critical mass for the PKS system: For M < 8ir, one has global solutions 
for which diffusion dominates so that all of the mass tends to infinity as the time tends to infinity, see [7] , while for 
M > 8n, solutions develop singularities, see [25] . 

Our focus in this paper is on the case M — 8ir. Notice that for M = 8tt, J-pks is exactly the functional that is 
on the left hand side in (|1.4j) . Since the densities q\ are minimizers of J-pks f° r M — 87r, it follows that 

S:Fp f [0x] (x)=O, 
dp 

and then from (|1.2p that each q\ - and each of their translates - is a stationary solution of (|1.2I) ; i.e., of 
Of course, this can also be checked directly. Our main goal in this paper is to determine the stability of these 
solutions, and to determine basins of attraction for them. In achieving this goal, we develop several novel functional 
inequalities, and a strategy of concentration control that may be useful elsewhere, and may be the main contribution 
of the paper. 

Each of the Q\ has an infinite second moment, and so shall all of the functions in the basins of attraction that 
we find for them. This must be the case according to previous work [6] on the case M = 8 7r for initial data with a 
finite second moment. The paper [5J proves the global existence of weak solutions with finite second moment that 
satisfy the free energy dissipation inequality 



f 




h 


Js. 2 



y ST PKS [p(t)) 2 
Sp 



d.r 



dt < Jpks[p(0)] , 



which is what one would guess should hold from (|1.3[) . Moreover, [B] proves that every such solution blows up at 
infinite time. That is, the t — > oo limit of any such solution is a Dirac mass 8 it at the center of mass of the initial 
data. Furthermore, a point mass of mass M is a stationary measure- valued solution in the sense introduced in |22j 
if and only if M > 8ir. Let us finally mention that an analysis of basins of attraction of g\ in the radial case was 
done in [4] in terms of certain relative moment conditions. 

From this point of view, the solutions in the critical mass case with finite initial second moment are choosing 
in their large time asymptotics the only possible stationary state with a finite second moment. 

Let us finally comment that the family of stationary solutions g\ play a role too in the conjectured profile of 
blow-up for any point singularity of the solutions for masses M > 8ir. Velazquez has proved (42j [43] that the inner 
part of the matched- asymptotics expansion for the blow-up profile is given by these stationary solutions for the 
critical mass value. 



1.2 The second Lyapunov functional 

The essential tool in our construction and analysis of solutions of the critical mass PKS system is an interesting 
and somewhat surprising interplay between the PKS system and another evolution equation which also has the q\ 
as stationary solutions - the Fokker-Planck version of the fast diffusion equation in R 2 with exponent 1/2: 



du I — 7T 

-£{t, x) = Ay/u(t, x) + 2 J &iv(x u(t, x)) i>0,i€l 2 , 
at V ™i (1 7) 

u{0,x) = u (x) > iel 2 , 

corresponding to the fast diffusion equation — = A\/u by a self-similar change of variable, see [41] . In the interest 
of brevity we refer to (|1.7|) as the fast-diffusion equation. 
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This equation can also be written in a form analogous to (ll.2[) : for A > 0, define the functional H\ on the 
non- negative functions in L 1 (M 2 ) by 



Hx[u] := f [VW)- VQ\(x)) 2 Qx~ 1/2 (x) dx. 



This functional is the relative entropy of the fast diffusion equation with respect to the stationary solution q\ . The 
unique minimizer of T-L\ is q\, and a simple formal computation yields 



6H\[u] _ _J_ J_ 
Su 

from which one sees that (|1.7[) can be rewritten as 

— (t,ar) = div f u(t,s)V — ^^0) 
It follows that for classical solutions u of (11.71) for which one can integrate by parts, 



(1.8) 



(1.9) 



dt 



n x [u(t)] 



u(t, x) 



,SH x [u(t)] 



(x) 



dx 



i V logu(t, x) + 2 Jjjj xy/u(t, x) 



dx 



As one sees from (|1.8|) and (|1.9|) , the densities q\ are stationary solutions of the fast diffusion equation (|1.7j) , as 
well as the PKS system (jl.lj) . This is much more than a coincidence, and there are very close connections between 
the two evolution equations. 

Unlike the functional J-pKSj the functional T~i\ is convex. In fact, it is not only convex in the usual sense, but is 
also displacement convex. This second type of convexity, arising in the theory of optimal mass transportation, will 
play a basic role in our analysis. We shall explain the relevant aspects of displacement convexity in Section 2 of 
this paper. For now, we return to the convexity of Hx in the usual sense. Making simple computations, one finds 
that 



H\[u] := y/nMX + J J^\x?u(x) -2y/v& + Vo>Jd 



dx . 



(1.10) 



Let us define the functionals 

Gi[u] 



\x\ 2 u{x) dx 



and 



02 M 



\J u(x) dx . 



Since Q\[u] is affine on its domain of definition, and since Q^lu] is convex on its domain of definition, one might 
formally conclude the convexity of Ti.\ on its domain of definition. In fact, those who are familiar with displacement 
convexity will recognize that functionals Q\[u] and Gi[v\ are displacement convex on their domains of definition. 

Unfortunately, separate consideration of Qi[u] and Q%[v\ is not relevant in our context: Note that \J Q\ix) is 
not integrable, and thus if "HaM is to be well defined, ■J u(x) cannot be integrable either. Furthermore, since 
|x| 2 gA(a;) is not integrable, it is clear that |x| 2 w(x) also will not be integrable on the whole domain of definition of 
"Ha- Thus, cancelations are crucial to the definition of "Ha, and the integral in (|1.10[) cannot be split into a sum of 
three integrals to be analyzed separately. 

As far as the convexity (in the usual sense) of "Ha is concerned, it is easy to give a rigorous proof: Indeed, H\[u] 
can be written as 

Hx[u]:= ( [<S>{u{x))-$>{qx(x))-<$>'{qx){u{x)- qx(x))] dx (1.11) 

with &(s) = —2^/s, which is a convex function. However, displacement convexity is essential to our strategy, and 
even after we have properly introduced the notion of displacement convexity, we shall have to work much harder 
to prove that T-Lx is in fact strictly displacement convex on its domain of definition. 

The convexity properties of T-L\ are relevant to the analysis of the PKS system due to the perhaps surprising 
fact that %\ is monotone decreasing also along solutions of the critical mass PKS system (|1.1|) . and not only along 
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solutions of the fast diffusion equation (|1.7[) . This gives us a second Lyapunov function for the critical mass PKS 
system. To see why this should be so, we make a formal calculation that we shall revisit in full rigor later on: Let 
p be a sufficiently nice solution of the PKS system. Then 



df 



5Hx[p ] 
Sp 

p\7 



div p(t,x)\7 



Sp 



dx = — 



Sp 



Sp 



dx 



Vex \fp. 



V [log p — G * p] dx 



AM 



V [log p — G * p] dx 



(1.12) 



Integrating by parts once more on the term involving the Green's function, 

|Vp| 2 



/ VVp-V[logp-G* P ]dx = l [ ^jr+/ v^AG*p=i/ -51-/ P 3/2 dx 

JR2 2 J R 2 p 6 ' 2 J R 2 2 J R 2 p A / 2 J R 2 



Also, L 2 xp ■ V log p dx — —2M and, making the same symmetrization that led to (11.6 



p(x) x ■ VG * p(x) dx 



-in 



x — y 

p(t, x) (x - y) ■ — — ^ pfa y) dx d v 



\x - y\ 2 



M 2 

47T 



(1.13) 



Using the last three calculations in (|1.12l) . we find 

4-Kx[pit)] = -\l l ^Sdx+[ p^dx + A 



Mtt 



8tt / 



Notice that the constant term vanishes in critical mass case M = 87r. Thus, in the critical mass case, formal 
calculation yields that for all T > 0, 



Kx[p(T)] 



|Vp| s 



R 2 P 



3/2 



(t,x)dx- / p 3/2 (t,x)dx 



dt < Hx[po] 



(1.14) 



In fact, the formal computation yields equality instead of merely inequality in (|1.14j) . but it is this inequality that 
is useful to us, and this is what we shall actually prove for the solutions that we construct here. 

The key to exploiting (|1.14j) is a particular case of the Gagliardo-Nirenberg-Sobolev (GNS) inequalities for which 
the sharp form was found by Del Pino and Dolbeault 20] . 

1.2 LEMMA (Gagliardo-Nirenberg-Sobolev inequality). For all functions f in R 2 with a square integrable distri- 
butional gradient V/, 



|/| b dx< / |V/| 2 dx / |/| 4 dz , 

JR 2 JR 2 

and there is equality if and only if f is a multiple of a translate of gx 1 ^ 4 for some A > 0. 
To apply this, note that at least for strictly positive densities p, 

|Vp| 2 



■(x) dx 



Therefore, we define: 



1.3 DEFINITION (Entropy dissipation functional). For any density p of total mass 8tt such that p 1 ^ 4 has a 
square-integrable distributional gradient, we define the entropy dissipation functional T>[p] by 

V[p}=8 f \Vp 1/i {x)\ 2 dx~ ( p 3/2 (x)dx, 
it 2 Jm. 2 

and we define V[p] = oo in all other cases. 
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1.4 LEMMA (Dissipation of W\). For all densities p of mass M = 8tt, 

V[ P ] > , 

and moreover, there is equality if and only p is a translate of q\ for some A > 0. 

Proof: Let / = p 1 / 4 and note that J R2 f 4 (x) dx — 8n. Multiplying V[p] through by n, the claim follows directly 
from Lemma ll.2l I 

1.3 The main results on the PKS equation 

The formal result (|1.14p may now be written as 

l-T 



H x [p{T)]+ f V\p(t)]dt<U x [p ]. 
Jo 

Since 'Hx{p{T)] > 0, this suggests at the very least that 

lim I f V[p(t)]dt< lim l<H x [ Po ]=Q 



and then Lemma 11.41 suggests that for all large t, p(t) must be close to for some p > 0. However, an easy 
calculation, see Remark 12.51 using the fact that J R2 \x\ 2 Q\(x) dx — oo, shows that 'H\{g^) — oo for /j / A. 
Therefore, since H\[p(t)] is non- increasing, one expects that p — A. In short, the formal calculations made so far 
suggest that for solutions p of the PKS system with initial data po satisfying H\[po] < oo, limt-j.oo p{x, t) = Q\(x). 
We now make one more definition, and then state our main results on the PKS equation: 

1.5 DEFINITION (Properly dissipative weak solutions of the PKS equation). Let Let po be any density on R 2 
with mass 8-7T, such that for some A > 0, H\[po] < oo. Let p : [0,oo) — > L 1 (R 2 ) satisfy: 

(1.5.1) For each t > 0, p(t) is a continuous curve of densities of mass 87r in the sense that for each bounded and 
globally Lipschitz function ip on R 2 , t >— > / ip(x) p(t,x) dx is continuous with p(0) = pq. 

(1.5.2) For each T > S > 0, and each smooth and compactly supported function ip on R 2 , 
ip(x)p(T, x) dx ~ J ip(x)p(S, x) dx + J J Aip(x)p(t, x) dx dt 

-if/ ^, P-. V " H -" , M ^.. 

4ttJs Jr2 xR2 \x-y\ 2 

(1.5.3) For each T > 0, 

U X \p{T)]+ [ V[p(t)} dt < Hx[p ] , (1.15) 
Jo 

so that p satisfies the entropy-entropy dissipation inequality expected for solutions of the PKS equation. 
Then p is a properly dissipative weak solution of the PKS equation (jl.ll) with initial data pq. 

1.6 THEOREM (Existence and regularity of properly dissipative weak solutions). Let po be any density on R 2 
with mass 87r, such that -^pksIpo] < °o> and for some A > 0, Hx[po] < oo. Then there exists a properly dissipative 
solution of the PKS equation with initial data po- Moreover, the solutions we construct have additional 
regularity properties, including: 

(1.6.1) For any S > and any p with 1 < p < oo, there is a constant C depending only on S, p, A and Hx[po] such 
that for all t > S, \\p(t)\\ p < C. 
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(1.6.2) The distributional gradient of p 1 ^ is square integrable over [0, oo) x R 2 . and in fact, 

[ [ \Vp^ 4 {t,x)\ 2 dxdt <Hx{p } ■ 
Jo Js. 2 

(1.6.3) J-pKs[p(i)] well defined for each t, and is monotone decreasing: J 7 pks[p(^)] < J~pks[p( s )] for all < s < t. 

1.7 THEOREM (Basins of attraction). Let po be any density on R 2 with mass 8ir, such that J-pksIpo] < oo, an d 
for some A > 0, rl\[po] < oo. Let p be any properly dissipative weak solution of the PKS equation (jl.l[) with initial 
data po satisfying the additional regularity properties (1.6.1), (1.6.1) and (1.6.3) of Theorem llM Then 

lim Tpks [/»(*)] = -^PKsteA] and lim \\p(t) - g\\\i = . 

t— too t— foo 

Let B\ denote the set of densities p Q of mass 87T for which Jpks[p] < °o, and H\[p] < oo. According to 
Theorem 1 1.71 B\ is a basin of attraction for g\ under the PKS evolution in the sense that any properly dissipative 
weak solution with initial data in B\, and the regularity produced here, converges strongly to g\ in L 1 (M 2 ). 

1.4 Controlled concentration inequalities 

The proof of the additional regularity in Theorem 1 1.61 and then Theorem 1 1.7[ might at first appear to be possible 
by a standard application of entropy-entropy dissipation methods, given the entropy-entropy dissipation inequality 
(I1.15|) . However, this is not the case. The essential point is that T>[p] is not a convex function of p, and even worse, 
it is a difference of two functionals of p that can each be arbitrarily large even when T>[p] is very close to zero. 
Indeed, for M = 8ir and each A > 0, T)[g\] = while 

lim || gA || 3/2 — oo , lim ||V£iA 1/4 ||2 = oo , and lim g x — 8irS . 

A— >0 A— vO A— >0 

the point mass of 87r at 0. It follows that the level sets of T> cannot be weakly compact in L : (R 2 ). Likewise, 
^pks[pa] = 87r(log8 — 1) for all A > while lim,\_j.o Qx — 8tt8q. Thus, a family of densities of mass 87r on which 
^pks is uniformly bounded need not be weakly compact in L 1 (R 2 ) 

In these examples of non-compactness for level sets of T> and .Fpks, we have a family of densities, which, in 
the limit, concentrate all of their mass at single point. We shall show here that this is essentially the only way 
compactness can fail for a family of densities of mass 8n on which T> or .Fpks is uniformly bounded. 

Compactness of level sets of T~L\ fails for more mundane reasons: A glance at (II. lip is enough to see that for all 
h > 0, one can construct a sequence of {pk} small perturbations of g\ lying in {p : H\[p\ < h} such that for some 
e > depending on h, liminffc^oo J,, pk{x) dx > e. Thus, level sets of H\ are not uniformly integrable, and 
thus not even weakly compact in L 1 (IR 2 ). 

However, the densities in level sets of H\ do have a crucial property: They must have "thick tails" and for 
this reason, they cannot concentrate more than a limited fraction of their mass on any given small set. Thus, the 
examples of non-compact subsequences in level sets of T> and .Fpks that we exhibited above do not lie in any level 
set of T-l\, and as we shall show, neither do any other non-compact subsequences. In Section 3 we prove: 

1.8 THEOREM (Thick Tails). Let p be any density of mass M such that H\[p] < oo. Then for ^ := \e~ 1/5 
and any s > 1 

/ p(x)dx>r)*e-^m nM [ g x {x) dx = ^h-e'vTm Ux[p] . 

J\x\ 2 >\s 2 J\x\ 2 >\s 2 1 + S 

Though the statement of this theorem makes no reference to optimal mass transportation, the proof we give 
in Section 3 relies heavily on the optimal mass transportation results we present in Section 2, including a new 
Talagrand type inequality, Theorem 12.41 involving %\. 

Using the very mild control on concentration provided by the Thick Tails Theorem, which says that densities 
p with H\ finite for some A cannot possibly concentrate most of their mass near any one point, we prove two 
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"compactness via controlled concentration" theorems. Of course, some sort of concentration control is inherent in 
any compactness theorem for densities, but the point here is that given only the very mild limit on concentration 
provided by the Thick Tails Theorem, we are able, in Section 3, to prove compactness for the level sets of J-pks 
and T>. The first of these theorems concerns J-pks- 

1.9 THEOREM (Concentration control for J-~pks)- Let p be any density with mass M — 8n, with Ti.x[p] < oo for 
some A > 0. Then there exist positive computable constants 71 and Cccf depending only on X and T-L\[p] such that 



7i 



/ p log + pdx < J"pks [p] + Cccf ■ 

JR 2 



Our second 'compactness via controlled concentration" theorem concerns T>: 

1.10 THEOREM (Concentration control for T>). Let p be any density with mass 8n, J-pks[p] finite, and T-L\[p] 
finite for some A > 0. Then there exist positive computable constants 72 and Cccd depending only on X, 7l\[p] and 
F~pks[p] such that 

72 / IVp 1 / 4 ! 2 dx < ttV[p] + Cccd ■ 
Jr 2 

Moreover, the same conclusion holds if we replace the assumption that J-pks [p] finite with the assumption that the 
entropy £[p] is finite, except that now the constant depends on the bound on £[p] instead of the bound on •FpksIp]- 

Theorems 11.91 and 11.101 give us the "vertical control" needed for a compactness result. The horizontal control 
is provided by Hx alone. Not only does a bound on T-L\[p\ ensure that p has thick tails, it also ensures that the 
tails are not too thick: A bound on H\[p] provides a bound on all moments of p up to but not including order 2. 
Unlike the Thick Tails Theorem, this result is elementary: 

1.11 LEMMA (Localization). For all densities p with mass M and all A > 0, and all < q < 2, there is an 
explicitly computable constant C depending only on q, X and M so that 



\x\ip(x)dx<C(l + H x {p}) q/2 . 

Proof: By the Cauchy-Schwarz inequality 

/ \p- £>a|£>a 1/4 dx < / \y/p- y/g\\gl 1/4 \y/p + y/gx\dx 
Jr 2 Jr 2 



< VH\[p]\\VP + V^h < y/2MHx\p] ■ 

Since g 3 ^ 4 is integrable, there is a constant C, depending only on A and M, whose explicit form is easily worked 
out, for which 

g- 1/4 pdx< [ Ql 1/4 gxdx+ x /2MHx[p} <C(l+Hx[p}) 1/2 ■ 



Now repeat the Cauchy-Schwarz inequality, but this time with p A 3 ^ 8 in place of g x then there is a (different) 
constant C, depending only on A and M such that: 

/ \p~ £>a|£>a 3/8 dx < / \^/p- y/g^\g^ 1/4 \y/p+ ^/g^\g^ 1/8 dx 
Jr 2 Jr 2 

< VnxJp\\\VpQx 1/s + ^2x 1/8 h < Vu^lp}xjc(i + Hxlp}) 1/2 . 

Since g 5 ^ S is integrable, then by changing the constant C accordingly, whose explicit form is easily worked out, we 
deduce 

gl 3/8 pdx< [ g- x 3/8 gxdx + C{\+Ux[p\f /4 . 
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The obvious iteration of this argument leads to 

/ Qx r P^< I er£Acb + C(l + HA[p]) 2r 

JR 2 JR 2 

for each r of the form 1/2 — (l/2) fc for k £ N, and then by interpolation, for all < r < 1/2, where of course C 
depends on r as well as A and M. To conclude, note that Q^ r (x) ~ |a;| 4r for large \x\. ■ 

Lemma \l . 1 II shows in particular that when H\[p] < oo, then plog(e + \x\ 2 ) £ L 1 (M 2 ), so that the Newtonian 
potential of p is well defined. Also, Lemma fl . 1 II together with Theorem 11.91 shows . via the Dunford-Pettis Theorem 
that the intersections of level sets of T-L\ and J-pks are at least weakly compact in L 1 (M d ), and stronger conclusions 
follows for the intersections of level sets of W\, J-pxs and T>. 

The "compactness via controlled concentration" provided by H\ and its dissipation T> through Theorem II. 101 
and Lemma 11.111 is the core of our proof of Theorem 11.71 However this is not the only use we shall make of 
compactness via controlled concentration: It is absolutely essential to our construction of properly dissipative weak 
solutions. 

Indeed, in many problems in which one seeks to prove an entropy-entropy-dissipation inequality such as ()1.15[) , 
both the entropy functional % and its dissipation T> would be weakly lower semicontinuous, often due to some 
convexity property. Then, if {p n }neN is a sequence of nice or approximate solutions of the evolution equation 
converging weakly to a weak solution p, one would have 



H[p{T)} < lim H[p n (T)} 



and / V[p(t)]dt< lim / T>[p n (t)]dt, 
Jo n ^°° Jq 



which is very helpful if one is trying to prove something like (|1.15jl . 

While in our case %\ is convex and lower semicontinuous, T> is the difference of two non-comparable convex 
functions and has no lower semicontinuity. Therefore, we need new tools to prove (1 1 . 1 5[) . and as we shall see, it is 
once again the compactness via controlled concentration that does the trick. 



1.5 Further developments 

One can build on the regularity results obtained here to prove additional regularity. Indeed, if p is one of the 
solutions we have constructed here, it is easy to prove that for any a > 0, Vc(x, t) is bounded and continuous on 
(a, oo) x R 2 , only using the continuity properties on p in t, the uniform control on first moments, and the fact that 
p(t) is uniformly bounded in both L 1 and L 3 for all t > a. Thus "freezing" b := Vc, p is seen to be a weak solution 
of the linear parabolic equation 

^ = Ap - div(6p) , 

with b bounded and continuous. Parabolic regularity theory may now be applied. In fact, the arguments developed 
in [TJ [2] can be applied to get L°°-bounds of the density for all positive times based on the L p -bounds, 
1 < p < oo, obtained in Theorem 11.61 A further development that requires new tools is to bound the rate of 
convergence to the equilibrium q\ in our convergence theorem. 

An interesting problem whose solution would lead to rate information is to characterize the stability of the 
GNS inequality that we have used. That is, we know that T>[p] = if and only if p is a translate of q\ for some 
A > 0, since, as we have seen, this is simply a restatement of a sharp GNS inequality of Del Pino and Dolbeault. 
A stability result for this inequality would be a result stating that, for any e > 0, if T>[p] is sufficiently small, then 
the distance, in some metric, from p to some translate of some q\ 1 A > 0, is no more than e. It would also be useful 
to quantify the qualitative stability result for the Log-HLS inequality that we prove and use in Section 5. Work in 
this direction is underway. 
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1.6 Other equations with a second Lyapunov functional 

The second Lyapunov functional Ti.\ is more useful to us than the primary Lyapunov functional J-pKSi which 
actually drives the evolution, because of its convexity properties, especially its displacement convexity, as explained 
in Section 3. 

There is a "canonical way" to produce gradient flow evolution equations that have a convex second Lyapunov 
functional that has been investigated in [32]. Indeed, both the PKS equation and the fast diffusion equation are 
gradient flow systems where the gradient is computed using the 2-Wasserstein metric, as we recall in Section 3. To 
keep things simple here, let us explain the mechanism studied in [32j in the finite dimensional Euclidean case. 

Let V be a smooth convex function on 1™. Let W be the smooth function on R" defined by W(x) — \ W(x)\ 2 . 
Now consider the evolution equation 

±x(t) = -VW(x(t)) . 

Then of course, for any solution x, 

±W(x(t)) = -\VW(x(t))\ 2 <0, 

and so W is monotone decreasing along the evolution. It is the primary Lyapunov function for this flow. Next, 
note that since VW = 2{RessV} W, 

— V(x(t)) = - [W • VW] (x(t\) = -2 { W • [HessF] W} (x(t)) < , 
at 

since the Hessian of V is positive. Thus, V is a second Lyapunov function for the gradient flow driven by W . 

An example in |32j concerns a porous medium equation on the line, which is gradient flow in the 2-Wasserstein 
metric for a certain entropy functional. With this entropy functional playing the role of V , the gradient flow 
equation for the functional corresponding to W is a certain fourth order equation of thin-film type. 

The fact that the entropy for the porous medium equations is a second Lyapunov functional for this fourth 
order thin film equation had been discovered earlier in |18) and exploited as the key to understanding the long 
time behavior of the latter equation. Again in this case, the second Lyapunov function is strictly and uniformly 
displacement convex, while the primary Lyapunov functional is not displacement convex at all. 

In the case studied here, the second Lyapunov functional does not arise through the mechanism studied in [31] , 
or any other evident natural mechanism, and we have no "explanation" of why one should expect H\ to decrease 
along the PKS flow. However, as explained in [32] . once one knows this, it is a consequence, formally at least, that 
Fpks decreases along the fast diffusion flow. This has interesting consequences that are investigated in [lOj . 

The motivation for doing the computation to check the monotonicity is twofold: First, both evolution equations 
have the same steady states, which is certainly necessary, but not at all sufficient, for the computation to work out. 
Second, there are many sharp inequalities that have negative powers of 1 + \x\ 2 as their cases of equality, so there 
are tools available to try to prove the positivity of the dissipation. 

1.7 A brief outline of the rest of the paper 

The rest of the paper is organized as follows. Section 2 begins with a brief summary of some results concerning the 
2-Wasserstein metric and gradient flows with respect to it. In particular, we recall a discrete variational scheme 
due to Jordan, Kinderlehrer and Otto j2S] for constructing solutions of a class of equations including both the PKS 
equation and the critical fast diffusion equation. We also recall McCann's [34] notion of displacement convexity, and 
explain how this should, at least formally, lead to the entropy-entropy dissipation inequality that we seek. Making 
the formal calculation rigorous will then be reduced to rigorously proving certain consequences of displacement 
convexity for T-L\, and this will be facilitated by the "robustness" of displacement convexity. 

The latter half of Section 2 is more novel. As we have noted earlier, T-L\ is formally a sum of displacement 
convex terms, however, for the densities that concern us, each of the terms is divergent. Thus, we are forced to 
introduce a regularization of While there are many tools available to regularize functions that are convex in 
the usual sense (e.g. infimal convolution), there is no general approach to regularizing functionals while preserving, 
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or at least not severely damaging, their formal displacement convexity properties. The regularization developed in 
the second half of Section 2 is one of the cornerstones of the paper. 

In Section 3 we prove the controlled concentration results that have been stated and discussed in previous 
subsections. 

In Section 4, we lay the ground work for the proof of Theorem 11.61 on the existence of properly dissipative 
weak solutions. These will be constructed using a variant of the Jordan, Kinderlehrer and Otto [26] scheme, which 
constructs the evolutions by solving a sequence of variational problems, as in di Giorgi's "minimizing steps" method. 

In this method, the Euler-Lagrange equation for the variational problem solved at each step often provides 
essential a-priori regularity on the minimizing density p. Once again, at this point in our problem, we encounter 
difficulties due to potential cancelation of infinities. To resolve these, we are forced to regularize .Fpks- The discrete 
scheme provides a very convenient framework in which to impose and control the regularization: We use a different 
degree of regularization at each discrete time step. Because of the regularization, we will at least know that at each 
time step, Vy/p is square integrable, but we shall have no useful quantitative bound on ||V*/jo||2' Still, this gives us 
enough regularity to make some crucial integrations by parts, and then eventually through the use of Theorem ll.101 
we shall obtain a useful quantitative bound on HVp 1 / 4 )^. 

In Section 5, we pass to the continuous time limit, and provide the proofs of Theorems 11.61 and 11.71 Here, the 
flexibility of choosing the degree of regularization at each time step is crucial to cope with the errors committed 
in the displacement convexity of the regularized functional and to get the right dissipation in the limit. The 
convergence in L 1 needs a qualitative control of the error in the log-HLS inequality: if J-pksIp] is small enough, the 
densities are closed to some in L 1 norm. This together with the fact that T-L\[p] is non increasing proves that 
p = A. 

2 Displacement convexity and the PKS system 

2.1 Gradient flows in the Wasserstein metric and displacement convexity 

We recall some facts concerning the 2- Wasserstein metric that will be used here. We shall be brief, aiming mainly 
to establish terminology and notation. For more background, see [44j[2]. Let V(M. 2 ) denote the set of probability 
measures in R 2 , and let V2 (R 2 ) the subset of probability measures with finite second moments. Define the functional 
W 2 in 7>(R 2 ) x V(R 2 ) by 

W 2 ( M , v) = inf ff \x- y\ 2 dH(x, y) , (2.1) 

lLel J JR 2 xB 2 

where II runs over the set T of all couplings of the probability measures p and v\ that is, the set of joint probability 
measures in R 2 x R 2 with first marginal p and second marginal v. For absolutely continuous probability measures 
/ da: and g dx we will simply write W 2 (/, g) in place of W 2 (/ dx, g dx). Clearly, W 2 is finite in 7> 2 (R 2 ) x 7> 2 (R 2 ), 
though it takes on the value +00 in certain pairs (jj,,v) G ^(R 2 ) x "P(R 2 ) - for example if p belongs to ^(R 2 ), 
but v does not. It is easy to see that W 2 is a metric on T-^R 2 ); it is called the 2- Wasserstein metric, where the 2 
refers to the exponent 2 on the distance \x — y\. More generally, given any v G T-^R 2 ), W 2 is a metric on the subset 
of -P(R 2 ) given by {p G V{R 2 ) : W 2 ((J.,v) < 00 }. 

A result of Brenier [8] as extended by McCann 33 , provides effective control over the minimization problem 
defining W 2 (/i, v). To recall this result, let T be a measurable map R 2 — > R 2 . We say that T transports p onto v, 
if for any measurable set B C R 2 , v(B) = p o T^ 1 (B). In this case we say that v is the push-forward of p by T, 
v = T#p. An equivalent formulation is that v = T^p if 

/ C(T(x)) dp(x) = f C(y) My) VC G C°(R 2 ) . (2.2) 
Jr 2 Jr 2 

By the Brenier-McCann Theorem [HI [33] , for any two probability measures p, and v on R 2 not charging Hausdorff 
dimension 1 sets, there is an essentially unique convex function ip in R 2 such that V(^#/i = v and 

W 2 (p,u)= f \x-V<p(x)\ 2 dp(x) . (2.3) 
Jr 2 
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The essential uniqueness is that if tp and tp are two such convex functions, then V^s = V</5 almost everywhere with 
respect to p. In this paper we will be concerned with densities whose mass is not necessarily one. If p and v 
are two positive measures of mass M > 0, we define W2(/i, v) in terms of the 2-Wasserstein distance between the 
probability measures p/M and v/M as follows: 

W!(m, v) = MW1(m/M, i//M) . (2.4) 

This normalization convention has the advantage that if Vip#(p/M) = {v/M), then (j2.3[) is still valid for arbitrary 
M. Note that if (|2.2p holds for /i and z/, it also holds if we change p and v by multiplying them by a positive 
constant, i.e., Vc^#(p/M) = {v/M) if and only if V(/?#p = 

In Section 5 we shall also use the p-Wasserstein distance, 1 < p < 2, especially for p = 1, on account of a 
useful description of compact sets for this metric. For two probability measures p and v on ]R 2 , p-Wasserstein 
distance W p (p, i/) is defined by (12.11) where 2 is substituted by p. For two positive measures of mass M, we define 
W p {/j,,v) = y/M~W p {(i/M,v/M). This normalization is chosen taking into account (|2.4[) to extend the standard 
ordering relation for the W p -metrics on probability measures; that is, by Holder's inequality, we have for any 
1 < P < 2 

W p (ji,v) <W 2 (p,v) . (2.5) 

A fundamental insight of Otto [35] is that the 2-Wasserstein metric is useful when considering any evolution 
equation on densities p that can be written in the form 



^(t,s) = div (/>(/.. r)V 



dp 



(2.6) 



for some functional Q. The prime example of (|2.6[) considered in [26] is the Fokkcr-Planck equation for probability 
densities for which 

Q\p\ = [ P {x) logp{x) dx + \ I \x\ 2 p{x) dx . (2.7) 

JR 2 z JR 2 

In [36] rescaled porous medium equations were also included. Otto's insight [35] [36] is that the equation (|2.6p is 
gradient flow for the functional Q with respect to the 2- Wasserstein metric. This is true for a large class of equations 
of the form flU), see [2j [Tl] [13l [16l EZ1 [44]. The "gradient flow in the 2-Wasserstein metric" point of view is 
useful to us here for two reasons: 



It provides the means for constructing well-behaved solutions of the equation in question through the solution 
of a sequence of variation problems; the Jordan-Kinderlehrer-Otto (JKO) scheme. 

It provides the means for studying the rate at which solutions of (|2.6|) converge to minimizers of the functional 
Q, at least when the functional Q has a certain convexity property. 

The convexity property referred to in the second point is McCann's notion of displacement convexity [34j . which 
turns out to be convexity in the "Riemannian metric" associated to the 2-Wasserstein metric; see [33]. If the 
functional Q is uniformly displacement convex, then there are automatically a family of functional inequalities that 
govern the convergence of solutions of (|2.6j) to minimizers of Q. In concrete terms, the functional Q is said to be 
displacement convex in case the following is true: For any two densities po and pi of the same mass M, let (p be 
the essentially unique convex function such that V</3#po = p\. For < t < 1, define 

\x\ 2 

ip t {x) = (1 - t)—* - + ttp{x) and p t = V</? t #p • 

The displacement interpolation between po and pi is the path of densities t <— > pt, < t < 1. Let 7 be any real 
number. To say that Q is ^/-displacement convex means that for all such densities po and pi, and all < t < 1, 
(1 — t)Q[po] + tQ[pi] — G[pt] ^ 7^(1 ~ *)Wa(po) Pi)- S is simply displacement convex if this is true for 7 = 0, and Q 
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is uniformly displacement convex is this is true for some 7 > 0. Let us recall the characterization of displacement 
convexity given by McCann in [34 i for functionals of the form 

Q*\p\ := / $(p(x)) dx, (2.8) 

where $ : [0, 00) — > [0,oo) with $(0) = 0. Then McCann's Theorem says that if s 1 — >■ s d $(s~ d ) is convex non- 
increasing on (0, +00) then the functional is displacement convex, and this condition is essentially necessary. 

A much simpler result, also from [34], is that if V is any real valued function on R 2 such that for all Xq, x\ € M 2 
and all < t < 1, (1 - t)V{x ) + tV(x{) - V({\ - t)x + txt) > jt(l - t)\x Q - xi| 2 , then the functional 

V[p] = [ V(x) p(x) dx 
Jr 2 

is 7-displacement convex. Using these results, one readily checks that in the case of the Fokkcr-Planck equation, the 
functional (|2.7| is indeed 2-displacement convex. The consequent inequalities that govern the long time behavior 
of solutions are Gross's logarithmic Sobolev inequality and the Talagrand inequality for Gaussian measures. Our 
analysis of long time behavior for the PKS system falls outside the scope of previous work in this direction since the 
functional J-pks is n ot displacement convex. The key reason that it is useful to bring the second formal Lyapunov 
functional %\ into the analysis of the PKS system is that it is displacement convex. In the next section we prove 
the displacement convexity of T-L\, and study its consequences. 



2.2 The critical fast diffusion equation as gradient flow of a uniformly displacement 
convex entropy 

The equation 

Ou 

— (t,x) = Au m (t,x) + Kdiv(xu(t,x)) , (2.9) 
at 

where k is a non-negative constant and to > 0, is called the porous medium equation with k = and m > 1 while 
for k = and < to < 1 is called the fast diffusion equation. When n > 0, there is a restoring drift. In case 
m = 1, (|2.9p is of course the heat equation for k = 0, and the linear Fokker-Planck equation for k > 0. 
Equation (I2.9P can be written in the gradient flow form 



^(t,x) = div ( u (t,x)V^-) with g\u] = I ( -^—u m {x) + K^-u{x) I dx , (2.10) 
at \ ou/ J R 2 \ m — 1 2 



which shows that the evolution equation (|2.9p is gradient flow for Q with respect to the Wasserstein metric. 

The value to = 1/2 for this equation in M 2 is critical in the sense that the functional Q in (|2.10[) is strictly 
displacement convex for to > 1/2, but is not displacement convex for to < 1/2. There are many of other "critical" 
values of to between and 1/2 at which other things happen, see [TH] for instance. But since displacement convexity 
plays a crucial role in our work, it is natural to refer to the to = 1/2 case as critical here. Indeed, by the criteria 
of McCann introduced above since Q\u] = V[u] + Q$ m [u] with V(x) — k\x\ 2 /2 and <p m {s) — s m /(to — 1), then V[u] 
is uniformly displacement convex, and for to > 1/2, t/^ m [u] is displacement convex. As might be expected, some 
difficulties arise at the critical value to = 1/2. 

Since (|2.9p is gradient flow for Q, one might hope to find stable steady states by finding the minimizers u of Q . 
Computing the Euler-Lagrange equation we find to/ (to — l)u m_1 + k|x| 2 /2 = C, where C is a Lagrange multiplier 
for the constraint M :— L 2 u{x) dx, which is conserved. In the case to = 1/2 and choosing 



K = KM > x:=2 \llTx wcfind u{x) = lF (x + \x\^ =gx{xh 

One readily checks that u — g\ is a steady state solution to (|2.9p with k = km. a and so the family of stationary 
solutions of the PKS system which we are investigating are also stationary solutions of the critical fast diffusion 
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equation for different drifts k = km,x- However, as neither ^J~q\(x) nor |x| 2 ^a is integrable, these functions are not 
in the domain of definition of Q , and so are not minimizcrs of 



Q[u} = J (^-2\J u{x) + Km,x^-u(x)J dx , 



(2.11) 



the m = 1/2 version of (|2.9p with k — hm,x as above. 

The cure is a simple renormalization as introduced in |191 130] : Consider instead the functional Hx [u] defined by 



u i V 



[y/u(x) - \Jgx{x)j dx + KM,x I -y- [u(x) - Q\(x)] dx. (2.12) 



Then, at least as long as u has the same behavior at infinity as does q\ 1 the integrals will converge. The counter 
terms that we have subtracted off from our functional do not depend on u, and hence they do not affect SQ/Su. 
This is the key idea used in the improvements of rates of convergence for the fast diffusion equation, see []j|J |3U] • 
Since km,aM 2 = ^1 \J~Qx ~ ^ k m.x, the functional in (I2.12[) can be written in the following simpler form, which we 
take to be the definition of the critical fast diffusion entropy: 



> 2 



H\[u}:= / - — '—dx. 

Vox 



It is easy to check that for in = 1/2, (|2.9[) can be written in the general form ()2.6[) with Q = Hx- As noted above 
the displacement convexity of Tlx is formally obvious from the fact that where y/u(x), \J Qx(x) and |a:| 2 u(a;) are 
integrable, H\[u] would differ from the right hand side of (12. lip by a constant. We provide a rigorous proof in the 
next subsection. 

2.3 Regularization of the critical fast diffusion entropy 

To show that u > Hx[u] is displacement convex, and more generally, to make rigorous computations involving 
critical fast diffusion entropy, T-L\[u], we introduce a regularized version of the critical fast diffusion entropy: 

2.1 DEFINITION (Regularized fast-diffusion relative entropy functional). For 5 > 0, and u a density with mass 
M, define Hx,s[u] by 



T-Lx.shA = \ 



(\fu + 5 - y/gx + S) ^ 



VOX + S 

2.2 PROPOSITION (Displacement convexity of relative entropy functionals) . For any density u € L^(M. 2 ) of 
mass M, 5 i— > Hx.s[u] is monotone increasing as S decreases to zero, and 

HmHx,jM =Hx[u] . (2.13) 

Furthermore, let Uq and u\ belong to L^R 2 ) of total mass M such that W2(ito,iti) < oo, and let u t , < t < 1 be 
their displacement interpolation. Then for each 5 > 0, 

(1 - t)H\,s[uo] + tHx,s[ui] -Hx,s[ut] > t(l-t)K s (u ,ui) , 

where Ks(uo,u\) satisfies 

lim K s (u ,ui) = k m .x W 2 (w ,wi) , (2-14) 

d— 

and K$(uq,ui) > j$ W|(uo,ui) with 75 < 0. Consequently, the maps u 1— > T-L\,s[u] are 75- displacement convex and 
the map u 1— > "Ha[m] is strictly uniformly displacement convex: 

{l-t)Hx [«o] + t H x - Hx [ut] > k m ,x t (1 - t) W 2 (u , u x ) . 
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Proof: As 8 n- ( \fu~+8 — y/gx + S) 2 is non-increasing. Then, as S decreases, 



Uu + 5 - + . (Vu- Jq\ 

increases to 



2 



V QX + S y/0\ 

By the monotone convergence theorem and (|2.13|) . the monotonicity in 5 follows. Next, 



- Vpa + 6) _ u _^ (f — I r ^ ( 5 



(VuTd - Vs) + - + + 5 -2V6 

v ' \ ox + o 



Vb\ + 5 Vpa + 5 V ' V~Q 

Where by the mean value theorem 

> <£ x 

- Vs 



< —= , + 8 — VS < — -= and 



Vq\ + s Vs' 2Vs 

These three terms are integrable and 



Vqx + S - 2 VS 



Vox+^ 

UxAu] = (I) + (II) + const. (2.15) 



where 

(I) 



/ — , u dx and (II):=2/ (VS~Vu + S) da 



The criterion (I2.8[) are easily checked for $(u) = VS — \/u + 8, and thus (II) is displacement convex. 

The term (I) in (j2 . 15[) is unfortunately not displacement convex in general. In fact, we will show that is 75- 
displacement convex with a explicit computable constant. In order to check the 7,5-displacement convexity of the 
regularized functional, notice that (gx + S)^ 1 ^ 2 is a function of \x\ 2 . Thus, the functional (I) is of the general form: 

V$(x) u(x) Ax 

with V$(x) = (gx + S)~ 1 / 2 . The characterization of 7,5-displacement convexity [5] ensures that this is implied by 
D 2 Vs > 7,5/2- Let us compute the hessian of the potential Vs(x). Define the function fs on [0,oo) by 



/ 5 (|x| 2 ) := Vs(x) - _ 
x/gx{x) + S 

We compute 

3AS(X + s) 



/*(*) = - 77 777U7- and f'M 



[A + S(X + s ) 2 f 2 JdVJ [A + 6(\ + s)*]W ' 

with A = M X/tt. Therefore 

D 2 V s (x) = 2f s (\x\ 2 )S l]+ 4f s '(\x\ 2 )(x®x), 
and taking into account that f'J (s) < 0, then 

£ • D 2 V s (x) ■ e > [2/KM 2 ) + Afl{\x\*)\x\ 2 ] \e := F s (\x\ 2 )\e (2.16) 

for all x, £ € IR 2 , where the function Fs is given by 

2A 2 + 2ASX 2 - 8AXSs - WASs 2 



F s (s)=2f s (s) + 4sf s '(s) = 



[A + S(X + s) 2 } 5 / 2 



It is obvious that the function Fs converges point- wise to the constant km, a as S — > in [0, 00). Moreover, since for 
each 8 > 0, the function Fs (s) — > as s — ¥ 00 and it is clear that is negative for s large enough since the denominator 
is positive and the numerator has a negative dominant term, then Fs attains its maximum and minimum in [0, 00). 
Then, we can choose its minimum value as 7,5 < and the 75-displacement convexity is proved. 
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In order, to show the limiting uniform displacement convexity, we need to refine our arguments. For that, we 
come back to the definition of convexity. Let tp be the essentially unique convex function such that V?/>#uo = u\. 
For < t < 1, define 



Vs(t) := I V s (x) [(1 - t)u (x) +tux(x) - u t (x)] Ax 

[(1 - t) V s (x) + t V s (Vip(x)) - V s (x + t(Vip(x) - x))] u (x) dx . 



We seek a lower bound on 7]g of the form r]s(t) > t (1 — t) Kg(u ,ui). Since 775(0) = r]g(l) = 0, it suffices for this 
purpose to show that r/g(t) > 2 Ks(u , ui) for all < t < 1. By denoting y :— Vi/)(x) — x, we have 



rfiit) = y- D'V 5 (x + ty) ■ y 1 u (x) dx . 

JR 2 

Using (|2.16l) . we readily obtain that 

»#(*)> / F s (\x + ty\ 2 )\y\ 2 u {x) dx > K s (u , m) , 

with 



K s (u ,ui) := min / F s (\x + ty\ )\y\ u (x) dx. 

0<t<lJ R 2 

Now, let us observe that the function F$ is bounded in [0, 00) uniformly in 8. For that, note that f' s is decreasing 
and thus f' s (s) < f' s (0) < A -1 / 2 . On the other hand, by the geometric-arithmetic mean inequality, we get 

M>^^ s " 2 aTJ?^' <217) 

As a consequence, we get 

g 

I I-F.5 Hi- (0,00) < ;^=|, 



and thus, 



\F s (\x + ty\ 2 )\\y\ 2 u Q (x) < —\ViP(x)\ 2 u a (x) G L\M. 2 ) , 



for all < t < 1. Thus, the dominated convergence theorem guarantees that 

Jim / F 5 (\x + ty\ 2 )\y\ 2 u (x) dx = K MtX \y\ 2 u {x) dx = k m ,aW|(u , Ui) 



8^0 



uniformly in < t < 1, which together with the definition of Kgiuo, u±) implies the uniform displacement convexity 
of the limiting functional T-L\[u\. I 

Continuing with the notation of Proposition 12. 2[ define the function hg on [0,1] by hs(t) = W\ y s[u t ] — 
Ks(uq, u\)t 2 . Then by Proposition 12.21 (1 — t)hs(Q) + ths(l) — hs(t) > 0, so that h$ is convex. Therefore, 
for alUe (0,1), 

Mi)-M0)> M "; M0) . 

This in turn implies that 

u r 1 v r 1 r ^a>M -UxAM , v , ^ 
ti\.8\u\\ - H.\.s\uq\ > hmsup h A5(u ,Wi) . 

t^o t 

To compute the limsup of the right hand side, we treat the two non-constant terms (I) and (II) in (|2.15[) 
separately. As we have noted (II) is displacement convex, and by well known theorems on the sub-gradients of 
displacement convex functions [U Chapter 10], this part contributes 

^~ "' ' (yip(x) — x) Uq{x) dx , 



2{u {x)+5f/ 2 
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as long as the integrand satisfies mild regularity properties; in particular whenever uo is bounded below on every 
compact set by some strictly positive number, and ^/uo has a square integrable distributional gradient. We shall 
show that both of these conditions hold in our application. Given that they do, then by the Cauchy-Schwarz 
inequality, 



/ 

JR 2 



|Vuo| 



(«o + S) 3 / 2 



\Vip(x) — x\ uq dx < 



IVkqI 
(«o + Sy- 



ria dx 



\Vip(x) — x\ 2 uo dx 

R 2 



|V^o| s 
uo 



dx W 2 (u ,ui) 



I [ |VV^| 2 dxW 2 {u ,u 1 ) 
Jk 2 



The contribution of (I) in (|2.15p can be treated by appealing to the general results in [5] since this functional 
is 75-displacement convex, in the notation of |2], and thus, this part contributes 

f' s {\x\ 2 )x-{V^{x)~x) u dx. 

Which is finite because s >-> y/s f$(s) is a rational functional which tends to when s goes to infinity. Thus we 
have: 

2.3 LEMMA (First-order characterization of displacement convexity). Let uo and ui be two densities of total 
mass M such that W 2 (uo,ui) < oo, and such that uq is uniformly bounded below on compact subsets ofM 2 by a 
strictly positive number, and that y/u~o~ has a square integrable distributional gradient. Let Vtp be the unique gradient 
of a convex function ip in M 2 so that V^#uo = u i ■ Then 



JR 2 



2Ax 



Vm 



(V-0(x) — x) uo dx + K$(uo,ui) 



\A+5{\+\x\ 2 ) 2 Y 2(uo + 5)i_ 
where Ks(u ,Ui) is defined in Proposition ^. 21 and the integrand in (|2.18p is integrable 
One might be tempted to take the limit 5 — > at this stage and to conclude 



(2.18) 



K\[uo] < U\[ui 



2u 



3/2 



(yip{x) - x) uo di-KA fi jW 2 (Mo,ui) , 



but without further information about Vip(x) — x, it is not possible to do this, or to justify the convergence of the 
integral. In our applications, it will be simpler to use the specific information that we obtain on \7?p(x) — x, then 
to do some integrations by parts, and then take the limit S — > 0. 

Let us finally deduce as an application of the uniform displacement convexity of the functional H\[u], an inter- 
esting functional inequality of Talagrand type. Actually, generalized Log-Sobolev-type inequalities lead formally 
to generalized Talagrand-type inequalities for this functional by repeating arguments due to Otto and Villani [381 
Theorem 1, Proposition 1] in the linear case and generalized in [16| Theorem 2.1]. Here, we are able to show it in 
full rigor by the previous approximation argument. 

2.4 THEOREM (Talagrand's inequality). Whenever u e £^(R 2 ) of mass M with H\[u] < oo, then 



km.x 



Proof: Using Lemma [2.31 with uo = g\ and u\ = u, we obtain that "Ha^M > Ks{g\,u). for all 8 > 0, since 

2Ax | Vg A _ Q 

[A + 5{\+\x\ 2 ff /2 2(Qx+S)V 2 
Thus, passing to the limit 5 — > taking into account (|2.14j) . the desired inequality is obtained. ■ 
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2.5 Remark (Basins of attraction). Theorem 12.41 tells us that the 2-Wasserstein distance of our initial data to 
the stationary state q\ is finite provided T-Lx[uq] < oo. Moreover, each of the equilibrium solutions Q\ are infinitely 



far apart in the W2 metric: We can easily check that with tp{x) 
uniqueness part of Brenier-McCann Theorem ensures 



'X/fi\x\ /2, one has V^p^Q^ — Q\- Thus, the 



1 



Qn{x) Ax = +00 



since the equilibrium densities Q\ all have infinite second moments. In particular, H\[g^] = +00 for p 7^ A. 



3 Proof of the concentration controlled inequalities 

3.1 Concentration control for J-pxs 

To prepare the way for the proof of Theorem 11.91 it is useful to give an elementary demonstration of a crude form 
of the log HLS inequality, without sharp constants, but which would nonetheless provide bounds on £[p] for all 
M < 8tt. 

3.1 LEMMA (Bounds on the entropy). Let p be a density of mass M on R 2 such that plogp and \x\p are in 
L 1 (]R 2 ). Then, for any a > 1/(8tt), there exists a constant C(M, a, A) > only depending on M, a and A such 
that 

~(G)+ * p(x) < a / p log p dx + C(M, a, A) + 4 a M log ( / + \x\ 2 pdx 
for al x £ M 2 . For a and X fixed, C(M, a, A) is monotone increasing in M . 

Proof: Recall the following Young type inequality: For all s,t > 0, st < slogs + e* _1 . Then, for any a > 0, we 
have 

st = a [s (t/a)] < as log s + ae t/a_1 . (3.1) 

We now apply this to 

* p(x) = -L / Qog\x-y\)-p(y)dy= f ^(- log \x - y\)-^rQx{y) dy , 
2 47rJ R 2 J\x- y \<i w Q\KV) 

under the integral sign with weight q\, and with 

■• — ' ' and t = — j- log \x — y\ . 



Q\{y) 4tt 
Since p\ is bounded above by M/(Xir), this yields 

The second integral on the right converges as long as a > 1/8tt, in which case, doing the integral explicitly, we find 

olM f 1 M 8na 2 

/ , tttta — 7 dz = for 8na > 1 . (3.3) 

ATrei^^i |^|V(^a) Ae87ra-1 v ' 

To relate the first integral to £ [p] , use the fact that s H> s log s is bounded below by — 1/e to conclude that 



P W^Ud y < / (^-)lo g (^-) ex dy + M 



\y-x\<l \Q\/ \QxJ Jr2 \Q\/ \Q\ 

< / plogpdy- / ploggxdy 



e 

M 

e 
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By Jensen's inequality for the concave function log in L 1 ((/9/M) da;), 

'AM \ 



/ plog q x dx = A/log ( — J - 4 / log (V^+ M 2 ) p da 

JR2 V f / JR 2 V 7 



> A/ log ^— ^ - 4M log 



Using ([33]) and (03 in (GO), we obtain 



1 ,ns i n ^ /" , i M 87ra 2 
-(G) + *p(x)<a p\ogpdy+—- - 

2 ,/ B 2 Ae »7ra — 1 



+ aM 



i - log T^J +31ogM + 41o g n ^/A + N 2 pda; 



The result follows with 



C(M,a,X) :=aM 



1 87ra 1 , / A , 

— + - - log - + 3 logM + 

Ae 8na — 1 e \ 7T 



(3.4) 



To apply this, let p be any density on R 2 of mass 8tt or less. For any number R > 1, define 

pi(x) := l { | 2; |>_ R} p(a;) and p 2 (a?) := l { | a; |<K}p(a;) . 

Also define 

pi(x) := l{| x |> fl _ 1} /3(x) and p 2 (a;) := 1{\ x \<r+i}P(x) ■ 
Then since the support of (G)+ has unit radius, 

/ p{G)+*pdx = I p(G) + *pidx+ p{G)+*p 2 &x 

JR 2 JR. 2 J R 2 

= / Pi (G)+ * pidx+ pi (G)+ * p 2 dx. 
Jr 2 Jr 2 

Now suppose it is possible to choose R > 1 so that for some < a < 8ir, 



/ pdx = 

J\x\>R-l JR 



pdx = / pi da; < 8ir — a 



and 



pdx — P2 dx < 8n — a . 



(3.5) 



'\x\>R-l JR 2 J\x\<R+l 

Then choosing a — (8tt — a/2) _1 , and applying the pointwise bounds from Lemma 13. II in (13.5[) . we obtain that 

1 f 167T 2d f 

-/ p(G)+*pda;< — plogp dx + 2C(8tt, a, A) + 32tt a log 

2 J R 2 167T - a J R 2 

It follows from this and the obvious fact that 



(3.6) 



\/\+Jx^pdx 



that 



— / pG*pdx>— p(G) + *pdx, 
Jr 2 ' Jr 2 

Fpks[p}>t7t^ / plogpda;-2G(8^,a,A)-32™logf / ^/X+Jxfpdx] . 

lbTr - a J R 2 \J H 2 J 



(3.7) 



If the mass M is less that 87r, the estimate p.7[) gives us an upper bound on the entropy of p in terms of ^pksIp] 
and the first moment of p. This would suffice, in place of the sharp logarithmic Sobolev inequality, to prove that 
no blow-up occurs or mass less than 8tt as it was done in [7] . This cutting- in-pieces argument is similar but simpler 
than the one in [6j Lemma 3.1] used for the bounded second initial moment case. 

If the mass is equal to 8n, one needs additional information to find an < a < 8tt and an R > 1 for which 
(|3.6p is true. The additional information in case the second initial moment is bounded was given by contradicting 
the convergence to a Delta Dirac, see [6j Lemma 3.1]. Here, we need to localize the mass and quantify the tails by 
using a bound on H\ [p] . One of the bounds is easy: 
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3.2 LEMMA (Solid core). Let p be a density of mass M such that for some X > 0, T-L\[p] < oo. Then 

/ , Pdx<— . (3.8) 

Proof: We start by reminding the bound 

\x\pdx< [ y/X + \x\ 2 pdx<2 y/\M + 2M 3 / 4 (A/vr) 1 / 4 ^% x [p\ , 

Jr 2 

proved in Lemma ll. Ill here written with explicit constants. The result is a direct consequence of this bound since: 

if / pdx<M/2, then / \x\p dx > rM/2 . 

J\x\<r JR. 2 

Choosing r — 4VA + 4:(X/Mn) 1 / 4 y / 'H\[p], we get a contradiction unless (|3.8[) is satisfied. ■ 
Now define 

and then using this value of R define p\ , p2 , pi and f>2 as above. By Lemma 13.21 

M , 
Pi ax < — < 47r . 

Thus, the left inequality in (|3.6I) is valid for this choice of R, and any a < ai :— 4ar. The inequality on the right in 
(|3.6p requires more work to achieve. The key is the following result on "thick tails" : 

3.3 THEOREM (Thick tails). Let p be a density of mass M such that Wf (p, q\) < oo. Then for 77* := |e _1/5 
and any s > 1, 

[ p(x)dx>r]*e-^ w *^ [ Qx (x) dx = ^ e -« ff ^«») (3.9) 

J\x\ 2 >\s 2 J\x\ 2 >Xs 2 l + S 

and 



[ p(x) dx > v „e~v^t nM [ 

J\x\ 2 >\s 2 J\x 



Bx {x)dx =^e-v£r^. (3.10) 
1 + s z 



(3.11) 



/|a;| 2 >As 2 J\x\ 2 >\s 2 

Proof: Let < 77* < 1/5 to be fixed later. Fix any s > 1 and define 77 by 

= I M 2> XS 2 p{x) dx 
V ' I\x\ 2 >xs 2 0A(aO da; ' 

We may assume that 77 < 77*, or else there is nothing to prove for this s. 

Let Vip be the gradient of a convex function such that V(p#Q\ = p. Define r := VXs, and define 

A := { x : \x\ > r and |Vy>(x)| < r } . 

Then since on A, |V<p(x) — x\ 2 > (\x\ — r) 2 , 

Wj(p,Q X )= [ \V^{x)~x\ 2 g x dx> [ (\x\ -r) 2 g x dx . 

JR 2 J A 

We now claim that 

r g x dx > (I-77) / gx{x)dx . (3.12) 

J\x\>r 
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Assuming (|3.12[) for the moment, 



W 2 2 (p, q x ) > mi { l(\x\-r) 2 g x dx : A C {\x\ > r} , Qx dx > {i - V ) / Qx (x) dx \ . (3.13) 

[JA JA J\x\>r J 

By the "bathtub principle", the infimum is achieved by choosing A to be the level set of (|a;| — r) 2 in {|a;| > r} that 
has the minimal mass, that is, A = {x : r < \x\ < r } where r is such that 

Qxdx = (l-r]) g\(x)dx, 

r<\x\<r J |:c|>r 

which is equivalent to 

Qx{x)dx = i] Q X {x)dx . 

\x\>r J \x\~>r 

By direct computation, this means MX/ '(X + r 2 ) = r/MX/(X + r 2 ), and then since 77 < 1, this means that r 2 > r 2 jr\. 
In fact, since 77 < 1/5, r > 2r, and so the optimal set A in (|3 . 13[) contains {x : 2r < \x\ < r/^/rj}. Therefore, 
combining the last estimates with (|3.13p . 

Wl(fi,Q\)> f (M - r) 2 g x dx > - f \x\ 2 gxdx. 

J2r<\x\<r/y/fj 4 J2r<\x\<r/^/rj 

Now recalling that r = V~Xs, explicitly calculating the integral yields 



W 2 (p,f?A) > 



AM 



log 



Remembering that s > 1 and t] < \ 
1 + s 2 /rj s 2 /rj 1 



> 



— and 



1 + s 2 /r] 
1 + 4s 2 



1 

4- - 



1 



rjj (l + s 2 /r/)(l+4s 2 ) 



1 



rjj (1 + s 2 /?7)(l +4s 2 ) 



> 4-- 



> 1 

77 J 5(1 + 77) ~ 5 



1 + 4s 2 5s 2 5?7 
Therefore, fixing 77* := e _1 / 5 /5 < 1/5, we get 

77 > r^e^™" " 2 

Combining this with (|3.1ip yields (|3 .9|) . Thus, to prove (13.91) . it suffices to prove (13.12[) . By definition, 



qx dx = I Qxdx> gx dx - / qx dx 

1 J{|a:|>r}n{|Vv3(a;)|<r} J|a;|>r ./|W(x)|>r 



Qxdx- I p dy = / Qxdx-ij q x dx 

|x|>r «/ I 3/ I > T J \ x\>r J I x I > r 



and this proves (|3.12[) . Finally, (|3.10|) follows from (|3.9[) and the Talagrand-type inequality, Theorem 12.41 
Now note that, by the definition of R, R 2 /X > 1, and hence Theorem 13.31 implies that 



Then since 



f p{x)dx>r,.^^Tt' Hx[p] ( Q X {x) 

J\x\>R+l J\x\>R+l 

I gx(x) dx > 

J\x\>R+l 



dx . 



!\x\>l 

8ttX 

W+T) 



Thus, with R as above and 

a 2 := 87^77*6 



H 



x/M^X I 



< 8tt 
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the second inequality in (|3.6[) is also satisfied. 

Proof of Theorem ll.9t Most of the work is now done since Q3.6P is satisfied by choosing a — min(ai, a 2 ). We need 
a final argument to take care of the control for the negative contribution of the entropy in terms of the localization 
of the mass of the distribution known as Carleman-type estimate. 

3.4 LEMMA (Control on the negative part of the entropy). For any density p e L\(R 2 ), if the moment 
J R2 m(x) p(x) dx is bounded with e~ m ( x > € L X (R 2 ) and m : — > Kq , then 

p(x) log_ p{x) dx < I m{x)p(x) dx + - [ e ~ m(x ^ dx . 

Jm 2 ' e j R 2 

Proof: Let p := px{p<i} an d M — jg2 f>( x ) dx < J R2 p(x) dx = M . Then 

p(x) (\ogp(x) +m(x)) dx = / [U(x)\ogU{x)]pdx - MlogZ 
where U := p/fi, p(x) — e~ m ^/Z with Z — J R2 e^ m ^ dx. The Jensen inequality yields 

/ [U(x) \ogU{x)}fidx > ( [ U(x)(j,dx) log ( f U(x)ndx) = M logM 

and 



p(x) log_ p(x) dx = / p(x) \ogp(x) dx > MlogM — MlogZ — \ m(x) p(x) dx 

Jm. 2 

> / m(x) p(x) dx . 

e J vd 



To control the negative part of the entropy in ()3.7j) . use Lemma 13.41 with m(x) = yA + |x| 2 and Lemma ll. Ill 

j olog_pdx> — / j olog_ ( o>— I m(x)p(x) dx / e ~ m ( x ) ^ x 



16ir — a Jb2 ' jjj2 Jjj2 e 

> - 2V\M - 2M 3/4 (\/ir) 1/4 y/n^p\ ~ ~ [ e~ m[x) dx 



e JR 2 

This gives the final control on the positive part of the entropy from (|3.7p : 

JtksW >t^ / /olog + pdx-2C(87r,a,A)- 327ralog( / y/X + \xf pdx] 

lb7r - a J R 2 \J K 2 J 

- 2V\M - 2M s / 4 (X/w) 1 / 4 v^aR- - / e~ m(a;) dx . (3.14) 

e Jr2 

Finally, we choose 

a 

T 1 := 

lo7r — a 

where a is given just above and an explicit expression for Cccf follows from (|3.14[) . I 
3.2 Concentration control for V 

Proof of Theorem 11.101 While in the proof of Theorem II .91 we used a "horizontal splitting" of p, here we use a 
"vertical splitting": Let / := p 1 / 4 . For /3 > 0, define f fj := minj/,/3 1 / 4 } and hp := f - f . We have 

ttV[p\=8tt [ \\7f\ 2 dx-n [ f dx . (3.15) 

Jm 2 Jr 2 
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Defining Ap = {x : f(x) > /3 1 / 4 } = {x : p(x) > /?}, we get 

fdx=f f$dx+ [ (hp + /3 1 / 4 ) 6 dx = [ f$ dx-0 2 \A f3 \+ I (hp + p 1 ^) 6 dx . (3.16) 

R 2 JR 2 \A 3 J A p 

By the convexity of x H> x 6 , for any 77 <E (0, 1) 

?3/2 



(hp +/3 1 / 4 ) 6 dx < L^\Ap\ + ^-L- j A h% dx. (3.17) 

By the inequality /J < y/^f 4 , and plugging (|3.17[) and (|3.16l) into (13.151) . we obtain 

TrV[p]>8n f |V/| 2 dx- n f h 6 p dx-Sir 2 ^-*^ 2 (\ - l) \A p\ . 
Jr 2 (l-^) 5 ^ 2 \V ) 



(3.18) 



By the GNS inequality, Lemma [1721 applied to hp: 



TV 



hpdx>- \Vhf}\ 2 dx(-. hidx) . (3.19) 



(3.20) 



~ jt^W J R 2 13 -~j M 2 l (Sl u ^(T^y R2 '^' 

By definition of fp and hp, V fp — in the support of hp so that 

/ |V/| 2 dx = / |V/^| 2 dx + / \Vhp\ 2 dx. 

JR 2 JR 2 JR 2 

Using (|37l9l) and (pHO)) in (pUg]) . we obtain 

irV[p}>8ir [ \Vfp\ 2 dx + f 8tt - 1 / hpdx)[ \\7hp\ 2 dx - 8ir 2 ^- ir (3 3 / 2 ( -1 - l) |^| . (3.21) 

We obtain a result of the type we seek under any conditions that ensure the second term on the right is positive. 
Our first approach uses Theorem II .91 and so requires that both .Fpks[p], and H\[p], for some A > 0, be finite. 
By (|3.ip once more and Theorem ll.91 we get for all a > 0, 



/ hidx= [ l Aff pdx<a( I plogpdx + e^-'lApl) <a ^ PKs[p] + CcCF +ae^ 1 \Ap\. (3.22) 

JR 2 JR 2 \JR 2 / 71 

This sort of estimate is frequently used in large deviations problems. Now choose 77 = 1/2 in (|3.21[) . and then 
choose a so that the first term in the right hand side in (|3.22l) is 7r/8, we have, for this choice of r\ and a, 



By Chebychev's inequality, | .A^l < and so we can choose (3 so that for any 7 G (0,47r), 

1 

(1 _ ^?) 5 



8?r ~ /-, -35 / h % dx ) - 7 ' 



as was to be shown with Cccd := 248/T 2 v / /3. Though we have explained how to compute /?, we shall not write 
down a formula. This proves the first part of the theorem. 

As for the second, note that we used the bound on J-pks only to obtain a bound on the entropy which was used 
in (|3.22[) . However, if we have by other means a bound on the entropy, we can use that in (|3.22[) in place of the 
bound on J-pKS- This proves the final part of the theorem. I 

4 Analysis of the discrete time variational scheme for the critical mass 
PKS system 



From now on, we will assume that the mass is 871- . 
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4.1 The Jordan-Kinderlehrer-Otto scheme for the critical mass PKS system 

The Jordan-Kinderlehrer-Otto (JKO) scheme for constructing solutions to (|2.6|) . as described in Section [21 would 
be to fix a time step r > 0, and inductively define the sequence {p k }ken by setting po to be the initial density, and 
then for k > 0, 

p^eargminj ^^ +gM} . (4.1) 

In other words, p k+l is some minimizer of the functional p M> W|(p, j o' c )/(2r) + Q[p\. Only existence of the 
minimizer is an issue, and not uniqueness, although in many examples that have been investigated a strict convexity 
argument furnishes the uniqueness. The key point is existence of a minimizer, since that provides a solution of the 
Euler-Lagrange equation for the minimization problem in (14.1| . Then, as shown in [26], the fact that each p k+1 
satisfies this Euler-Lagrange equation means that, under certain conditions on Q, if one defines p T by 

p T (t,x) = p k {x) for kr < t < (k + 1)t , with p T (0,x) — po , 

there is a sequence of values of r tending to zero along which p T tends to a solution of (|2.6[) in a suitable weak sense. 
This scheme of constructing weak solutions of the PKS system for M < 8ir was developed in [5]. However, for 
AI = 8tt we can not proceed in a very direct manner. Our problem lies outside the scope of previous applications 
of the JKO scheme, since at the critical mass M — 87r, (|1.5p provides no upper bound on £[p], and hence, it is not 
even clear that minimizers exist for the variational problem in (|4.ip when Q — J-pks an d M — Sir. Our controlled 
concentration inequalities could be used to solve this problem, but other more thorny issues arise when we would 
try to analyze the Euler-Lagrange equation. 

To circumvent these difficulties, we introduce a regularized functional. In fact, for reasons that will become 
evident later on, we shall even be forced to choose a different degree of regularization at each time step. 



4.2 Regularization of J-pks 
Let 



7 (x) := ^e-N 2 / 2 



be the standard Gaussian probability density in R 2 . Then, for all e > define %(x) = e~ 2 j(x/e), and define 
the regularized Green's function G c = 7 e * G * 7 e , where * denotes convolution, and G(x) = — l/(2n) log |x|. The 
radially symmetric, C°° probability density 7^ is the fundamental solution of the heat equation satisfying for 
t > 0, x e M 2 : 

d 1 A 

^7V7 = 2 A7V7 ■ 

It follows that 

d 1 1 d 

g- t G *^Vi= 2 AG * 7 ^ = " 2 7 ^* ' and thuS ' d~e Gt = ~ 2e72£ ' (4 ' 2) 
The right hand side is strictly negative everywhere. From this we deduce a useful monotonicity property: 

£l <e 2 => G ei (x) >G e2 (z) (4.3) 

for all i£l 2 . Let us point out that 

J^G £ = -2 72e - 8e 2 A 72e = 2x ■ V 72£ + 2 72c = 2 (l - l2e , (4.4) 
since j e satisfies e 2 A7 e + div (xj e ) = 0. 

4.1 LEMMA (First properties of G e ). Let G e be defined as above then: 
(i) For all x 6 M 2 ; G e (x) < G(x). 
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(ii) There exists C > such that for all x G R 2 , G e (x) < Ce' 2 . 
(Hi) For all (x,y) G M 4 , 

G t (x - y) > -i- [4+ log (e + |x| 2 ) + log (e + |y| 2 )] . 

Proof: (i) As 7 is radially symmetric and subharmonic in R 2 so that, by the mean value property the first item 
holds. 

(ii) Since log_ |x| is locally integrable in R 2 , for any i£l 2 



/ G(x - y)7 e (y) dy < ^- / log_ \x - y\ 7 £ (y) dy < 

J\y\<e Z7T J\y\<t 



G * J e (x) 

since 7 is bounded. Thus, we get 

G e (w)<^ f le (z) dz = - 2 . 
e Jv- e 2 

(iii) From the elementary inequality \z — w\ < \z\ + \w\ < 2max{|z| , \w\}, we obtain 

log \z — w\ < log 2 + log \z\ + log \w\ . 

Therefore, 

G(z -w)> --^(2 + log \z\ + log H) ■ 
Integrating both sides against 7 e (a; — z)j e (y — w), and using Jensen's inequality, we find 



G e ( X -y)>-- 



2 + log / \z\ lt (x-z)dz +lo. 



> 



2tt 



\w \ r y e (y — w) dw 



2 + log \z\j e (z)dz +log [\y\+ / |w| 7 e (iw) dw 



> t 4 + l0 S ( e + W") + l0 S ( e + l»| 2 )] . 

at least for e small enough so that J R2 \z\j e (z) dz is small enough. I 

One of the main uses that we will make of the regularization of the self interaction functional is that it provides a 
regularized density for the chemical attractant: given a mass density p, we define the regularized chemical attractant 
density c e by c e (x) = G e * p(x). 

4.2 LEMMA (Uniform estimate regularized chemoattractant). For all e > and all densities p with mass 8ir, 
the regularized chemical attractant density c e = G e * p satisfies 



and 



MVcelU < 8C HLS IMMIMtIU/s + 4+ ^ ||7l| 2 / 3 IINpII 



(4.5) 



Here Chls denotes the constant of the sharp Hardy-Littlewood-Sobolev (HLS) inequality [3T] for the special 
case p = q = 4/3: 



/(*)] 



1 



\x-y 



■g(y) dx dy < Chls II / 1 1 4/3 115 1 1 4/3 



(4.6) 



Though the explicit value of Chls is simple enough, see [3T] , our bounds and their proofs will perhaps be easier to 
read if leave Chls unevaluated in them, as a marker of the use of the HLS inequality. 
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Proof: By the Young inequality, we have ||Vce||oo = ||VG e * p\\oo < 8n ||VG e ||oo- And by the HLS inequality, 

|| VG.Iloo / 7e(- " Z) 7^-| le(w - y) dz dw < ^ \\ 7e \\l /3 = || 7 ||i /3 . 



Using the triangle inequality |x| < \x — z\ + \z — w\ + \w — y\ + \y\ we have 

||a;|Vc e | < 2 |(|x|7 e ) * VG * 7 e * p| + |7 e * p\ + \j e * VG * 7 e * |a;|p| . 
Using the Young and HLS inequalities, see (|4.6I) . we obtain 

GhLS 2 

IlklVCell^ < 8G H LS ||Fl7e|| 4 /3 IMI4/3 + 4 + "^T" IWI4/3 IIfHIi . 

and the second part of the result is obtained by using UTe II4/3 = £ ^ 1 ^ 2 |) T 1 1 4/3 - 

Using the regularized Green's function G e , we introduce the regularized self-interaction functional W e 

We(p) = // p(x) G e {x - y) p(y) dx dy . 



4.3 LEMMA (Continuity of the regularized interaction energy). Let p\ and P2 be any two densities in R 2 of mass 
8vr bounded in L 1 (M 2 ,log(e + |a;| 2 )). Then, for all < e < 1, 



\W e [pl]-We[p2}\ < 



- + 2Ce- 2 

7T 



||/3l||i 1 (R 2 Jog( e +|x| 2 )d.T) ||Pl ~ P2|Ul(BL 2 ,log(e+|x| 2 )dx) (4-7) 



Moreover, let {p n }neti be a sequence of densities all bounded in L 1 (M 2 ,log(e + \x\ 2 )) uniformly in n. If {p n }nevi 
converges weakly in L 1 (M 2 ) to p, then for each e > 0, 

We[p] < liminf We[pn] ■ 

n—>oo 

Proof: By definition of G £ : 

W e [p] = W [p * 7J = / |VG * p * 7e | 2 (z) dz > . 
Jr 2 

Using p :— p± — p2 we write 

W £ [pi] - W e [p 2 ] = W £ [p] + 2 / Pl * G e (ar) p(a:) dx > 2 f [p* G e ] pdx . (4.8) 

JR 2 JR 2 

Then combining Lemma 14.11 and (|4.8j) , we obtain 



We[pi] - We [P2\ > 2 / / p{y) G £ (i - y) p(x) dydx + 2 p{y)G e {x-y) p{x) dy dx 

J J p>0 J Jp<0 

+ 2 ~2 II p(y)p( x ) d v dx 

e JJp<0 



> - 



3 „G 
- + 2- 

7T 



||Pl II i 1 (M 2 ,log(e+|^| 2 ) da) IIpII L 1 (K 2 ,log(e+|x| 2 ) dx) 



Now swapping the roles of p\ and pi, we obtain (|4.7p . 

By Lemma 14.11 {p^}„ew bounded in L 1 (R 2 ,log(e + |a;| 2 ) dx) uniformly in n implies that p n * G £ is bounded in 
L°°(R 2 ) uniformly in n. Since {p„}„ e N converges to p weakly in L 1 (R 2 ), then p n * (xr G e ) — > p*(xn G e ) point-wise 
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for given any cut-off function \R with support in 5(0, i?) and thus p* G e € L°°(IR 2 ). Therefore, applying 
with p\ =: p and p2 '■= p n , we have 

liminf (We [p n ] ~ W e [p\) > Km 2 f [p* G t ] (p n - p) dx = , 

where we have used the weak convergence on the right hand side. I 
We are now ready to introduce our regularized free energy functional. 

4.4 DEFINITION (Regularized free energy functional). For all < e < 1, define 

•TpKsM : = / p(x)logp(x)dx- - [ p(x)G e (x~y)p(y) dxdy 

Jr 2 z JJr 2 xr 2 

on the set of densities p of mass 87r such that p G L 1 (M 2 , log(e + \x\ 2 ) dx) and plogp is integrable. 

The following lemma is an immediate consequence of (|4.3[) : 

4.5 LEMMA (Monotonicity in e). For all densities p of mass 8ir such that p <E L 1 (R 2 , log(e + \x\ 2 ) dx) and plog p 
is integrable, e i— >• -^pksW * s monotone decreasing in e. 

Note that by Lemma [4.11 (ii) and (iii), p(G e * p) is integrable for p s L 1 (M 2 ,log(e + \x\ 2 ) dx). Moreover, by 
Lemma 14.11 (i) 

^ksW > JtksW ■ (4.9) 
In particular, by the sharp log HLS inequality, see Lemma 1 1.1 1 

JpksW > -C(8tt) = 87r(-l + log8) . (4.10) 

By Lemma [33] (iii) , we have the upper bound independent of e: 

■TpKsM ^ / P( X )1°EP( X ) dx + 32?r + 2||p||ii(R2 ilog ( e+ | 

x\ 2 ) dx) • 

(4.11) 

4.6 LEMMA (Error estimate for regularized free energ y). For all p € L\ n L 3/2 (M 2 ) loit/i mass 87T, and a// 
e < (2Ve)-S 

JpksW ~ -Tries [p] < G 7 1 1 p| 1 3^ e . 
Proof: We use Holder's inequality and Young's inequality for convolutions to get 

■FPKSH-^PKSM <|Ml2 /3 ||G e -G|| 2 . 

Holder's inequality gives .Fp KS [p] — -^pksIp] < V8~7r II Pll 3/2 II — G|| 2 - Since G € L* oc (R 2 ), then G e (a;) — > G(x) a.e. 
in M 2 as e — )• 0. By Lemma [4.11 G c < Gg < G for < e < e. This implies that ||G e — C| 1 2 is non decreasing in e 
and it has zero limit as e — >■ by monotone convergence theorem. Moreover, the same arguments also show that 

lim||G e -G e -|| 2 = ||G e -G|| 2 . 
On the other hand, using (I4.2[) and (I4.4[) for |x| > 2e, we get for any fixed x € M 2 



G;(i)-G«(i)<(£-e) 



< 2e 2 7e (x) 



where second-order Taylor expansion at e of the function G e (x) for < e < e was used and the second order term 
is nonpositive due to (|4.4[) . Taking the limit e — > and integrating, we deduce that 

\ 1/2 

(G(x)-G e (x)) 2 dx <2e|| 7 || 2 . 

\x\>2e I 
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Since G £ < G, we can thus directly compute 

/ (G(x) - G t (x)) 2 dx < 4 f \G(z)\ 2 dz = Wire 2 (\ - log(2e) + (log(2e)) 2 ) . 

J\x\<2t J\z\<2t \ l J 

Finally, simple computations show that when | log(2e)| > 1/2, the term in parentheses on the right in no greater 
than 5 1 log(2e)| 2 . Collecting all together leads to the result with G 1 explicitly computable. ■ 

4.3 Existence and first properties of the JKO scheme minimizers 

Let S denote the set : 

S :={pe L^R 2 ) : / p{x) dx = M, W 2 (p, q\) < oo, £ [p] < oo, / |a;|p(x) da; < oo} . 
Jr 2 " Jr 2 

By dUm, the functional 

P ^ nZ + -^pksH 

is bounded from below on S. The next lemma asserts that it has minimizers, and begins the task of their analysis. 
We state this lemma for a single step since we shall be changing the value of e from step to step. 

4.7 THEOREM (Existence of minimizers). Let A > 0, < r < 1 and < e < 1. For all p e S 

argmiJ^%^l + ^ KS [p]l (4.12) 



pes { 2r 

is not empty, and each minimizer p belongs to S. Moreover, there exists K\ > depending only on £[po] and 
J R 2 \x\po dx such that Fp KS [p] — -hi- 
proof: Let {p^}keN be a minimizing sequence i.e. such that 

k—foo y zt J pes y It J 

By what we have noted just above, the infimum on the right hand side is finite. The following observation is the 
starting point for obtaining all of the bounds we need: Considering the trial function p = po itself, one sees that 
we may suppose 



- r pK S [P (fe) ] < 2V 2 7 PU; + ^PKSM = ^PKSIPO] 



Wl(pW,po) + ^ Kg[p(fc)] ^W|(po,p ) 
for all k. Consequently, for all k, 

-^PK S [P (fe) ] <-^pks[Po] and W 2 (p^,po)<2T[^ KS [p }-^ KS [p^}] . (4.13) 

We first bound W2(p^ k \ Q\) uniformly in k. Since po € S, (|4.11j) ensures that J- PKS [po] < oo, and provides a 
bound depending only on £[po] and J R2 \x\po dx. Then (|4.10[) provides a universal lower bound on FpusiP ]> an d 
thus by (|4.13p . there is a finite constant K\ depending only on £ [po]j and J R2 \x\po dx such that for all k, 

^ KS [p ik) ] < Kx and W 2 2 (pW,p ) < K x ■ (4.14) 

In particular, by the triangle inequality, for all k, ^N2{p^ k \ Qx) < VKx~ + W2(po, Qx) < oo. 

We next bound the first moments of p^ uniformly in k. Let V<p be the optimal transportation plan X/ipftp^ = 
g\. Then since |x| < \x— \7(p(x)\ + \ V<p(x)| for all x, integrating against p^ and using the Cauchy-Schwarz inequality 
yields 

/ |x| /><*>(*) dx< V8iW 3 (p (fc) ,eA)+ / \x\ex(x)dx< v / 8^[ V / ^ + W 2 (p ,^)] + / \x\g x (x)dx. 
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The right hand side is finite and independent of k. 

We next bound £[p^]. By part (ii) of Lemma T4.ll there is a constant C such that 

logpWO*) dz < ^ KS [p^} + < J| KS [po] + < Jfi + , 

6 6 6 

where we have used (|4.14[) once more. Again the right side is finite and independent of k. 

The last two uniform bounds show that {p^}keN is uniformly integrable. Hence, by the Dunford- Pettis theorem, 
there exist a weakly in L 1 convergent sub-sequence whose limit we shall denote by p. 

By a standard weak lower semicontinuity argument (see e.g. [2] for the weak lower semicontinuity of Wf ), p 
satisfies each of the three bound that we have proved uniformly for {p^}keN, and thus p e S. 

It remains to prove that the functional -^"p K g is lower semi-continuous on L 1 (K 2 ). For the entropy part, this is 
standard. For the self interaction part, this follows from Lemma |4~31 So that 



J£ KS [p] < liminf ^ KS [p (k) ] < K x . (4.15) 
Finally, the weak limit p is a minimizer. I 

4.8 PROPOSITION (Strict positivity of the minimizers). Let po satisfies the conditions of Theorem 14.71 Then 
any minimizer p of (|4. 1 2[) is uniformly bounded below on compact sets, i.e. for all R > 0, there exists s > such 
that 

p(x) > s almost everywhere in Dr :— {x : \x\ < R} . (4-16) 
Moreover, s does not depend on the chosen minimizer of Problem (|4.12j) in case of non-uniqueness. 

4.9 Remark (Idea of the proof). Let us try to quantify this simple statement: Entropy abhors a vacuum. The 
functional derivative of £[p\ is log p. On any set where p is very close to zero, logp is very negative, and we can 
lower £[p] by transporting some mass from where p is relatively large to this spot where it is very small. This will 
lower the entropy by a very large multiple of the transported mass. On the other hand, if we do not have to transport 
the mass too far, the effects on W 2 (p, po) and W e {p) will be relatively small. 

Proof: For any s > 0, let 

Ar(s) := {x £ Dr : p(x) > s} and Cr := / p(x) dx . 

Ja r (2/R2) 

For any s > 0, let ctR(s) :— {x £ Dr : p(x) < s}, and let |aij(s)| denote its Lebesgue measure. By Theorem 14.71 
|x|p(x) is integrable, and hence 



/ p(x) dx = / p{x) dx — I p(x)dx>8ir— j 
JD n Jk 2 J\x\>r Jm 2 



\x\ 

■ ! — L p(x) dx > 4 7T , 

R 



as long as AttR > L 2 M/°(aO dx. 

If |q;_r(s)| = for some s > 0, there is nothing to prove: p is bounded below uniformly by s on Dr. Therefore, 
suppose that |q;_r(s) | > for all s > 0. Pick some small positive numbers S and s, and define a new density p by 
transporting a mass S Cr |oir(s)| from Ar(2/ R 2 ) to o.r(s), distributing it uniformly there, which raises the density 
there by SCr. In formulas, choose s < 2/R 2 to have cifl(s) n Ar(2/R 2 ) = 0, and define a new density p by 

\l-6\a R (s)\)p(x) xeA R (2/R 2 ), 
P( x ) = { p{x) +5C R x£ a R (s) , 

p(x) otherwise . 

In order to ensure positivity, we have to impose 5|q;_r(s)| < Sir R 2 < 1/2. In this way, it is easy to check that p is 
a density. 
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Note that \\p — p\\ i < 2 8 \c(r(s)\Cr , and since all the modifications take place on Dr, 

||p-p||z,i(R^iog(e+H 2 )daO < log(e + R 2 ) 2 5 \a R (s) | C R . 

It now follows from the bounds on p derived Theorem 14.71 and from Lemma l4.3l that there is a constant K depending 
only on R, e, £[po] and H\[po] such that 



We[p\<We[p} + S\a R (s)\K . 
Using Taylor's expansion of x i— > xlogrr, that log a; is increasing and assuming s < 8Cr, we obtain 



(4.17) 



/ [p\ogp- p\ogp]Ax = / 
Jr 2 J a r 



(2/-R 2 



{(l-6\aR(s)\)plog[(l-6\aR(s)\)p]- plogp} dx 



[{p + 8 C R ) log (p + 8 C R ) - p log p] dx 



< S\a R ( S )\ 
<S\a R (s)\C R 



p log [(1 - 8 \a R (s)\) p]dx + S \a R (s)\ C R log(2<5 C R ) 



-log( — ) + \og{28C R ) 



(4.18) 



where 8 \o.r(s)\ < 1/2 and x G Ar(2/R 2 ) were used in the last estimate. 

To estimate the difference W^p, po) — W 2 (p, po), let II denote the optimal coupling of p and po, and use it to 
define a non-optimal coupling II of p and po. To do this, let p be the measure supported on Ar(2/R 2 ) with density 
p, and hence total mass Cr. Let v be the dx-uniform distribution on olr^s) with total mass Cr. Let \7tp be the 
optimal transportation plan with Vip#p = v 1 and define the map T : M 2 — > R 2 by 



T(x) = 



Vip(x) x G A R (2/R 2 ) 
x otherwise . 



Then II, given by II = (1 — 5|afl(s)|)II + <5|a fl (s)|(T (g> Id)#II is a coupling of p and p : an d hence 



W^(p,p ) < 



|x - y\ 2 dn(x,y) 

(R 2 

(1- 5\a R (s)\) W 2 (p,po)+S\a R ( s )\ 



\T(x)-y\ 2 dIl(x,y). 



Then, since \T(x) - y\ 2 < 2\T(x) - x\ 2 + 2\x - y\ 2 , and \T(x) — x\ < 2i?, since all of the transportation induced by 
V"0 takes place inside Dr, it follows that 

W 2 (p, P() ) < (l + 5\aR( S )\)W 2 (p,p ) + 8\aR(s)\(8ir) 2 8R 2 . 

By bounds on p derived in the proof of Theorem 14. 71 there is a constant K depending only on £[po] and Hxlpo] 
such that W 2 (p,p ) < Kt. Finally then, there is a constant depending only on R, r, £ [p ] and W\[po] such that 



W 2 2 (p, Po ) <W 2 2 (p,p )+5\a R ( S )\K 
Combining (|4T8|) . (|4~T71) and (|4T9|) yields 



(4.19) 



+ -^pksN < + •FpksIpJ + <5 |a fl (s)| Cfl 



1 



log ( ) + log(2JC fl ) + X' 



2r ™ 5Lrj " 2r 

with a given constant -FT'. If |or(s)| > for all s > 0, then choosing 6 small enough such that 



log 



1 

R 2 



+ \og{28C R ) + K' < 
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contradicts the optimality of p. For instance, choosing sr = S/Cr, the procedure described above can be carried 
out, and we conclude that p is bounded below by sr on Dr. This proves (|4. 16[) . ■ 

We now continue the analysis of the minimizers p begun in Theorem 14.71 We obtained p G S and the lower 
bound (|4.16[) directly from the variational principle, but to proceed, we need the Euler-Lagrange equation for the 
variational problem (|4. 12|) . 

By the Brenier-McCann Theorem, there is a a lower semi-continuous convex function ip in R 2 such that 'S/ipffp = 
po, and V<p is uniquely determined on the support of p, which is all ]R 2 by (|4.16|) . The Euler-Lagrange equation 
for ()4.12|) relates p, po and V</s: 

4.10 LEMMA (Euler-Lagrange equation). Let po satisfy the conditions of Theorem 14.71 and p be any minimizer 
for the variational problem in (|4.12|) . and let V<p be the unique gradient of a lower semi-continuous convex function 
such that Vipjfcp = po- Then the distributional gradient of p satisfies 

id-Vtp . 
-Vp + pVc e = -p (4.20) 

T 

where c e = G £ * p. In particular, since c e is differentiable everywhere, and ip is differ entiable almost everywhere, p 
is differentiable almost everywhere. 

The proof follows exactly the original procedure in [35], see also [5J Theorem 3.4], and we skip it here for the 
sake of brevity. The interested reader can see the details in the preprint version of this paper. 

4.11 LEMMA (Qualitative regularity estimates). Let po satisfy the conditions of Theorem 14.71 and let p be any 
minimizer for the variational problem in (14.121) . Then yfp has a square integrable distributional gradient, and for 
any 1 < p < oo, p p is integrable. 

Proof: By the positivity of p, see Proposition 14.71 we can divide both sides of (|4.20[) by y/p, to obtain 



where V</? is such that Vtp^p = po- By the triangle inequality, 



1/2 1 , , N 1/2 

- , ..; A„\ ill U Y7,„^M2 

1/2 



2||VVp|| 2 <( j^\Vc e (x)\ z p{x)dx) +-[ I _\x-V<p(x)\ J p(x)dx 



<[ I \Vc e (x)\ 2 p(x)dx) +-W 2 (p,p ). (4.21) 



R 2 



1 



T 



By Lemma 14.21 IIVc^oo is uniformly bounded, and so is the first term of (|4.21|) . This proves that ^fp has a square 
integrable distributional gradient. The integrability of pP is then a consequence of the following classical version of 
the GNS inequality valid for functions on R 2 with p € [2,oo) 



\v\ p dx < D p 



IVd 2 d£ 



p/2-l 

v\ 2 dx 



/K 2 

applied to v = ^J~p. I 

4.12 Remark. Although the bounds in Lemma 14.111 are not quantitative, and would certainly be e dependent if we 
were to extract quantitative bounds, we shall use them only to justify certain integrations by parts, and otherwise 
show that substraction of infinities does not invalidate computations that follow. Thus, these qualitative estimates 
are all we require concerning V ' J~p and p. However, they are absolutely crucial for their purpose, and their necessity 
is the main reason we have had to introduce the regularized Green's function G e , and along with it, the regularized 
chemical attractant. Without the regularization, we would only know that 2Vy / p— Vc^/p was square integrable - 
but the possible cancelation effects would not allow us to conclude that V \fp was square integrable. 



BCC August 1, 2011 



33 



4.4 A discrete form of the entropy-entropy dissipation inequality 



Our main goal in this subsection is to prove a discrete version of the entropy-entropy dissipation inequality (|1.15|) . 
The key idea is to use the k x displacement convexity of %\ and the "above the tangent" inequality for convex 
functions as follows: For given initial density po, let p be any minimizcr for the variational problem in (I4.12[) . Let 
Uti < t < 1 denote the displacement interpolation between p and po starting at p and ending at po- Then uq = p 
and u\ — pq. Since T~i\ is displacement convex, the "above the tangent" inequality for convex functions says that 



n x [ P ] + -n x [u t ] 



KxW 2 2 (p,p )<n x [p } 



A formal computation of the second term on the left would give, for e = 0, 

= V[p] . 



t=0 



Indeed, assuming Lemma 12.31 holds for S = e = applied to uq = p and u\ — po, we get 



H x [p] <Hx[p ] 



K X X 



Vp 



P 



3/2 



(V<p(x) - x) p dx - k x Wj 0, po) 



Using (|4.20[) . i-e. (V<p(x) — x)p — r (\7p — pVc) and expanding, we can rewrite this as 

2 



U x \p] <H x [po}-- 



ivpr 



2 p 



.3/2 



dx — 



Vc • Vp 



dx — k x / x -Vcpdx + k x / x ■ \7p dx 



k x Wf(p,p ) 



■=n x [p ] - ^ [(I) + (II) + (HI) + (IV)] - Kx W 2 (p,Po) . 
Using — Ac = p we have 

(II) = -2 / Vc • V dx = -2 / p 3/2 dx . 
Jr 2 Jr 2 
Using the symmctrization argument we obtain 

(III)=K A y-( / pdx) = 167TK A . 

And by integration by parts (IV) = — 16ttk x , resulting into 'Hx[p} < %a[po] — t^[p] — k a Wf(p, po)- 

However, to do the calculation in a rigorous manner we must take into account that e > 0, and we must use the 
regularized entropy functional "Ha,<5- Before proceding with this, we point out that no such estimate can be given 
for .Fpks since this functional is not displacement convex. 

4.13 LEMMA (Convexity estimates at the regularized level). Let po satisfy the conditions of Theorem 14.71 and 
let p be any minimizer for the variational problem in (|4.12[) . then 



H\,s[p] <"Ha,«5[po] - ^ 



|Vp|' 



dx + T 



(p + S) 3 / 2 

16ttk x t + 2C £ || \2f' s - k x \ (1 + |a;|)p||iT 
Ks(po,p) , 



p 3/2 dx + 16 7T v^r + 16nVS( J 7 + 1 + C* £ )r 



167T J,, 



2A 



r + 2r / V-[x/K|x| 2 )]p(x)dx 



and 



n x .s[p] <n x . s [po] - J / . dx + r / p 3 / 2 dx + 16 tt v^t + 16ttV5( J 7 + 1 + C e )r 

2 Jr 2 (p + d) 3 ^ J R 2 

- 16™ A r + 2C e || |2/i - « A | (1 + IxDpllxr + %/i || |x| tIU/3 \\p\U,*t 

V 2A 

+ 2r / V-[x/^(|x| 2 )]p(x)dx-/^(p ,p) , 

JR 2 
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where K$ is defined in Proposition 12.21 and the constants J 7 , C e; and C e are explicit constants. 

Proof: This is an elaborate calculation in which a number of integrations by parts operations must be carefully 
examined for boundary behavior. It is relegated to the Appendix. I 

As a consequence of this lemma, letting S go to 0, we obtain the following result concerning the dissipation of 
T-L\ in one discrete time step. 

4.14 COROLLARY (Convexity estimates). Let po satisfy the conditions of Theorem 14.71 If p is any minimizer 
for the variational problem in (|4.12[) then 

Hx\p]<Hx\po]-TV\p]+TA\\i\\ A/s -K X WZ(p,po) , (4.22) 

and 

Hx[p]<Hx[p ]-TV[p]+TV~eA\\p\\ 4/3 - Kx W 2 2 (p,pa) . (4.23) 
where A := 32tt(2A)- 1 /2 Chls 7H4/3. 
Proof: Let us first observe that 

\7 ■ [x f s (\x\ 2 )]p(x) dx = f [2f s (\x\ 2 ) + 2\xff s '(\x\ 2 )]p(x)dx. 

Js. 2 

Let us recall from the proof of Proposition [2~2l that 2f' s (s) /• n\ and f'g'(s) —> as i5 — > for all s > 0. Moreover, we 
have that sf$(s) is a bounded function uniformly in S from (|2.17p . These properties together with the dominated 
convergence theorem leads easily to 

II |2/J — kaI (1 + |a:|)p|U -> and / V • [x f' s (\x\ 2 )}p(x) dx -> 8k a tt 

Js, 2 

as S — > 0, since (1 + \x\)p G L 1 (M 2 ). By monotone convergence theorem, we obtain 

hm/ -^Ldx=t raCd*. 



8^0 J m2 (p + (5) 3 / 2 ' 7 R2 p 3 / 2 

Putting together all these facts and Proposition 12. 2[ we can pass to the limit as 6 — >• in Lemma H. 131 to get 
the desired estimates (|4~2"3"]) and P~2"2"J) . ■ 

4.5 One-step estimates 

Neither of the one step dissipation estimates that we have so far, namely (|4.23|) and (|4.22j) . are exactly what we 
need. The problem is the term 7"^4|| ^/|| 4/3 in the first of these, and the term t"^4.|| />| 1 4/3 in the second of these. 
These terms might be large compared to the other terms so that these estimates might even give only "negative 
dissipation" . 

In the first main result of this subsection, we use one and then the other of these inequalities in combination 
with the controlled concentration inequality of Theorem 11.101 to produce the kind of dissipation estimate that we 
really want. In the second main result, we show that LP norms of the densities are essentially propagated along 
each step of the discrete variational scheme. Again, Theorem 11.101 plavs a crucial role in both proofs. 

4.15 THEOREM (One-step theorem). Let p satisfy the conditions of Theorem 14.71 p be any minimizer for the 
variational problem in (14.121) . and choose any C po such that 

n x [p ] < C pa . (4.24) 

Define Q > 0, r„ > by 



Qq ■= C pa - H\[po] and Tq := min < , 1^ , (4.25) 
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where A is the constant given in Corollary 14.141 Finally, given Qq and < r < Tq , and also any positive integer 
t, let q be given by 



,1/3 



;^/3^ 7 -2/3 (7rCpo+r , CccD) 2/ 



Then for all r < Tq and all e = ei, p satisfies 



/'i i 



and 



Hx[p] - Hx[ P o] < - rV\p] + Q^t 2 2- 1 - « A W 2 (p , p) . 



Proof: By (|4.22l) , our choice of r and Qq in (|4.25[) implies that 

Ux\p] < n x [po] - rV[p] + 9±= Cpo -Q a - rV[p] + 9l< Cpo - rV[p] 



(4.26) 
(4.27) 

(4.28) 



On one hand, the GNS inequality, see Lemma [L2l implies T>[p] > so that (|4.28l) implies that p also satisfies (|4.24l) . 
On the other hand, since Jixip] cannot be negative it implies 

V\p\ < ^ . 

T 



Moreover, by the monotonicity of e .7-p KS [p], Lemma 1431 and (I4.15[) . J-pksL°] < -^pksW < +°°- We can thus 
apply the concentration controlled inequality, Theorem 11.101 which implies 



Vp 1 / 4 

By the GNS inequality of Lemma 11.21 we have 



2 dx < — [ttV[p] + Cccd] < - — [ttC Po + t*C C cd] 

72 T 72 



Vp 



1/4 



1 



c 3 



dx < [7rC Po + Tq C C cd] := — ■ 

r 72 t 



(4.29) 



T 



2/3 



p 3/2 dx < 8 

Next, by Holder's inequality, 

/ p 4 / 3 dx = [ p 1 ' 3 pdx < (Sn) 1 / 3 ( [ p 3 ' 2 dx) V3 < (Sir) 1 / 3 ( 

Jr 2 Jr 2 \Jr 2 J V 

Now using this bound in (|4.23[) . we obtain 

Ux[p] -Hx[p ] < -tV^+t 1 ' 3 ^ [a(8tt) 1/3 C 2 3 /3 ] - K X W 2 2 { Pa ,p). 

We thus obtain the stated result by choosing e — ee for any positive integer I. I 

4.16 LEMMA (Propagation of the L p -norm). Let po satisfy the conditions of Theorem \A.l\ Assume additionally 
that po €E L P (R 2 ), 2 < p < oo, and let p be any minimizer for the variational problem in (|4.12l) . then there exists 

Ko > which only depends on / p\ log p| dx such that for all K > Kq 

Jr 2 



[ (p-K) p + dx< [ (po~K) p + dx + TA 1 +TA 2 V[p\ 

Jr 2 Jr 2 



where A\ and A 2 are universal positive constants depending on K . 
Proof: The displacement convexity of the functional 



p^r [ ( P -K)\dx 

Jr 2 
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with 2 < p < oo and K > 0, is easy to check using McCann's criterion (I2.8[) . The Euler-Lagrange equation of 
the variational scheme (x — V<p) p = — rVp + rpVc e together with the standard first-order displacement convexity 
characterization [44J [2] imply 



(p - Kf, dx (p - dx < -p 



(p-A)^ 1 (V(,5-x)pdx 



<-(p-l)pr[ \V{p-K) + \ 2 {p-Kf- 2 dx 
Jr 2 

+ (p - l)r / V [(p - Vc e dx + prA / V [(p - Kf' 1 



Vc e da; 



. 4(p-l) 

< T 



{p-Kfl 2 dx + (p-l)r/ (p- A)^(-Ac e )dx 



prA / (p - Kf-^-Ac.) dx := r(h +h+ J 3 ) 

JR 2 



(4.30) 



The last two integration by parts have to be justified for any given e working as in the proof of Lemma 14.131 in the 
Appendix. Integrating by parts on the ball of radius R, we obtain for any k £ {p,p — 1} 



f V[(p-K) k + ]Vc e dx= [ 

J\x\<R J\a 



(P-K)l(- 



\x\<R 



< 



\x\<R 



Ac e )dx + (p - K) k + V 'c £ ■ n da 

J\x\=R 

f p k da. 

J\x\=R 



(p- K) k + (-Ac e )dx+\\Vc e \\ c 



It is enough to show by dominated convergence theorem that there exists a sequence of radii {iij }j e pj such that 
the boundary terms tend to zero as j — > oo. Due to Lemma 14.111 with p > 2, for any given natural N > 1, we can 
write that 



y / / p k da dr = / p k dx < oo , implying that lim / / p k da dr 

N=1 Jn-1 J\x\=r JR 2 N ->°° Jn-1 J\x\=r 



0. 



for k 6 {p,p ~ 1}, and the two integration by parts for any given e are justified. 

We now estimate Ii and ^3, showing in particular that they are finite. Starting with I2, using — Ac £ = p e where 
p e := 7 e * p * 7 £ , so that by Holder's inequality and Young's inequality for convolutions, obtain 

(4.31) 



(p - Kf + {-Ac t ) dx= (p - K)lp e dx < || (p - K) + \\l +1 \\ Pe \\ p+1 < \\(p A)+||£ +1 ||p|| p+1 
'■ Jr 2 

Likewise for I3, we use the fact that on the support of (p — A) + , A < p. Therefore 

K f (p- Kf^i-Ac,) dx = A f (p- K f-'p, dx < f (p- K) p - 1 Pe p dx . 
Jr 2 Jul 2 Jr 2 

Therefore, by Holder's inequality and Young's inequality for convolutions, 

A [ {p-K)\-\-Ac e )dx<\\{p-K) + \%l\\\pf p+1 . 
Jr 2 

Applying the arithmetic-geometric mean inequality to the right side of (|4.3ip . we have that for any v > 0, 



(4.32) 



Making a similar estimate for the right hand side of (|4.32p . and combining results, we have that 



h + h < FMWiP K) + \\lX\ + F 2 (i/)||p||£i 



11P+1 



(4.33) 



BCC August 1, 2011 



37 



where F\(v) is a positive linear combination of negative powers of v, and F2(v) is a positive linear combination of 
positive powers of v. 

By Lemma T4. Ill p p is integrable for any 1 < p < oo, and so the right side of (|4.33|) is finite. Then by (|4.30j) . 



{p-Kff 



dx < oo . 



From here we show that 

V (> 2 ) 



dx < 



Indeed, 



and 



V (> 2 ) 



dx 



(p-Kff 
(p-K)f 



2 1 Q K (2p-l)/2 

dx + [ttV\p] + Cccd] 



72 



(4.34) 



/ 




){P<K} 





dx < 16K (2p -^^ 



1/4 



dx 



dx < 



'{p<K} 
72 



V (V /2 ) dx , 
irV[p] + C CCD ] 



where we applied again the concentration controlled inequality, Theorem 11.101 using ()4.26|) . 
Following an idea of Jager and Luckhaus [25] . we use the GNS inequality 



v p+1 dx<D p ( I |Vv p/2 | 2 da; 



v dx 



(4.35) 



which is a consequence of the Sobolev embedding inequality ||u|| L 2( R 2) < C|| Vu||ii(R2) applied to u — i/p +1 )/ 2 and 
Cauchy-Schwarz inequality since Vu = p ^—v 1 / 2 Vv p / 2 . Applying (|4.35p to v = (p — K) + , we get 



(p - K) p + +1 dx < M(K) 
Then (|4~33l becomes 

h + h < Ft{v)M{K) 



{p-Kfl 2 dx where M{K) := I (p-K)+dx 



(p K) P J 2 



dx + F 2 (v) 8tt 



Vp 



p/2 



dx . 



We finally work with I\ to estimate it using (|4.34p as 



2(P-1) 



h < 



(p K)f 



dx 



Vp 



J'/ 2 



dx + 



mK (2 P -l)/2 
72 



[itV\p] + Cccd] 



Now choose vq > small enough such that 871-^2 (z^o) < 2(p — l)/p, and then Kq < 00 large enough such that 
M(K)F\(vq) < 2(p — l)/p. This choice of Kq only depends on vq and the bound on J R2 p\ logp| dx since 



We find 



M(K)= [ (p- K) + dx< [ pdx<—^— [ p logp dx < — !— [ p\og + pdx 
Jr 2 J P >K \ogKJ p>K logK J R 2 

/ (p - Kf + dx - (p - Kf + dx < r- ^ >- [irV[p] + C CC d] 

Jr 2 Jm 2 



for all K > Kq. The desired result follows with 

A\ = OcCD 

Pl2 



Pl2 



32if(2p-i)/2( p _ 1) 
and A 2 = — -n . 

PI2 
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5 Proof of the main results 
5.1 Approximate solutions 

We now combine the single step operations described in the previous section to inductively dchnc infinite sequences 
{(O^jfegN giving a discrete-time approximation to the PKS evolution For the rest of this section, fix any 

A > 0, and any density p in R 2 with total mass 87T with p log po integrable, and such that there exists C Po 
with H\[po] < C Po . It then follows from Lemma 1 1 . 1 1 1 that \x\po is integrable, and from the Talagrand inequality, 
Theorem 12 .4[ that W2(poj Qx) < °°- Thus, po satishes the conditions of Theorem 14 . 71 on the existence of minimizers 
for our single step variational problem. 

Fixing an arbitrarily small parameter r > 0, we now inductively dehne the sequence of densities {p^}kefi with 
p® := po by solving the sequence of variational problems 

p k T G argmin Py' 1 ) + 7& s \p]\ (5.1) 



pes y 2t 

for a sequence of regularization parameters {efejfcgN to be specified now. By Theorem 14. 71 the sequence {p k }keN is 
well defined no matter how we choose < r < 1 and {efcjfegN. We moreover define 



C P0 - HxlPr] > for each k . (5.2) 



5.1 LEMMA (Good step sizes). Let {p£} fceN be the sequence of minimizers defined inductively using (|5 . 1 1) starting 
from p T — Po- With Qk defined as in (|5.2p . let A be the constant given in Corollary 14. 141 and let A be defined by 

A: =n 0-V 

m=l v 

and note that 1 > A > 0. Choose any r > satisfying 



r<mhW ,1}:=t*, (r,.8) 



2A|| 7 ||4/3 



and define eu by 



,J/3 4-2/3 , n r , ^2/3 Qo _2 -* /v „\ 

.7T ' A7 2 (7rCp +CccD) J = — t I . (5.4) 



Then for all k, Qk > AQo > 0. In particular, 

AToie i/iai /or some constant Z, e/. '■= Zt 10 ^ 3 4~ fc , 

Proof: We shall show by induction that for each positive integer j 

j 

Qj " 



II 1 - -r ) Qo ' (5 - 5) 



which is somewhat more than we need since the right hand side is larger than AQo. 

We now make the inductive hypothesis that for some positive integer k, (15.51) is true for all positive integers 
j < k. Since A < 1, we may apply Theorem 14.151 with p^T 1 in place of po, and p k in place of p and with r and 
specified as above. Then the conclusion (|4.27[) can be simplified and rewritten as 



Ih-x r,-k 

n\[P T \ ^ n\[p T j t •/ - 

Since r < 1, this means that 



n,[p k T ]<n x [p^ 1 ]+r 2 ^2- k . (5.6) 



Qk > Qk-i ( 1 - — 



2J«- 

By the inductive hypothesis, we obtain (|5.5|) for j = fc. The proof that (|5.5|) is valid for j = 1 is a direct application 
of Theorem 14.151 in the same way, since A < 1 . ■ 
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5.2 The passage to continuous time 

Throughout the rest of this section, we assume that < r < t* , where t* is defined in (|5.3|) . and that €k is defined 
by (|5.4[) . and then that {p k }ken is a corresponding sequence of minimizers of (|5.ip . 

We now interpolate between the terms of the sequence {p k }keN to produce a function from [0, oo) to L 1 (M 2 ) 
that we shall show to be, for sufficiently small t, an approximate solution of the PKS system. For technical reasons, 
we shall need two distinct, but closely related, interpolations. 

• The Lipschitz interpolation: For each positive integer k, let Wip k be the optimal transportation plan with 
S7ip k #p k = p k ~ x - Then for (k — l)r < t < kr we define 

Mt)=(^-^Id + ^V/)#^. 

• The piecewise constant interpolation: For each t and each positive integer k with (fc — l)r < t < fcr we define 

p T (t) = Pr' 1 , With p T (0) = p . 

For displacement convex functionals of p, such as T-lx[p\, £[p], or the absolute first moment, any uniform bounds 
on the functional along the sequence {p^j/cgN extend to p(t) for all t, since if Q is such a functional, then for 
(k - 1)t < t < fcr, 

G[Pr{t)\ < tzStzlkgtfri] + k -I^±g[p k ] . 

T T 

Of course it is evident that for any sort of functional 0[p], displacement convex or not, a uniform bound on Q[p] 
along the sequence {p^}k^n extends to p(t) for all t. Some of the functionals with which we work, such as T>[p], 
are not displacement convex, and this is the reason we need the second interpolation. 

The uniform equicontinuity properties that we prove next explain the utility of the first interpolation, and also 
why we can use the two different interpolations at once. Since p\ is a minimizer for (|5.ip . using p^" 1 as trial 
function yields 

and hence, 

W 2 2 (p k T ,p k T -i) < 2r [J^s^-J^sM . (5-7) 
and by the monotonicity of e n- J-p KS [p] see Lemma 14.51 

7-pKs[^]<-FfWPr]<+°°- (5-8) 

In standard applications of the JKO scheme, in which the functional in the variational problem does not change 
from step to step, one would sum both sides in (|5 . T[) over a range of values of k, and then the sum of the terms 
on the right would telescope. This is not so in our case. However, for small e, J-p K g ~ -^pks an d we recover the 
telescoping sum in a useful approximate sense. The precise version of ^pks ~ -^pks follows from (|4.9[) . Lemma 14.61 
and (|5.4p . which says that e k = Zr 10 / 3 4:- k to get 

F&M- 1 ] -^pksIp^ 1 ] <C- ( Z\\p k -^ 2 r w /H- k < Z\\p k -'\\f 2 r*2- k (5.9) 

for t < f* := min(Y*, (2Z v / e) -3 / 10 ) with Z := C 7 Z according to Lemma WM We thus deduce 

W 2 2 (p k ,p k - 1 )<2r (^ PKS [p k - 1 ]~^ PKS [p k })+2Z\\p k -X / / lr i 2- k . 

Using (|5.8[) and (|4.29p as in the proof of Theorem 14. 151 where the concentration control inequality (| 1 . 10[) is crucial, 
we deduce that ||/3t _1 |I3/2 — C3/ 7 " to conclude that 

WK/^p*- 1 ) < 2r (-FpksIp^ 1 ] -^pks[^]) +2ZC 3 r 3 2~ k . (5.10) 
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We are almost in a position to obtain a crucial a-priori Holder continuity estimate, but there is still one more 
consequence of our step dependent regularization to deal with: If for each k we had been using the functional 
J-"pks instead of the functional -Fp^g, ^ would be immediate that k J-"pks[Pt] would be decreasing. Since by the 
Log-HLS inequality, J-pks is bounded below, this would give an immediate upper bound on the sum of the right 
hand side of (|5.10p over any range of k. 

However, we have used our freedom to choose the sequence {efcjfceN of regularization parameters to converge to 
zero as rapidly as we may require, and hence easily obtain: 

5.2 LEMMA (Uniform bounds on the free energy .Fpks)- There are positive constants J-q, T\ depending only on 
the initial data and the regularization mollifier 7 such that for each t < f* and each k G N, 

Proof: Directly from the variational problem (|5.1|) we have ^pksIPt] ^ -^pks^t Then, as above from (|4.9|) . 
(jOl) . and (Ei~2"9l we get 

•FpW^ 1 ] < ^ks iP^ 1 ] + ZC 3 r 2 2- k . (5.11) 
This means that the free energy Jpks is almost decreasing along {p*}ke®- A telescoping sum argument yields 
^PKstPr] — -^pksI' 00 ] + ZC 3 t 2 , an( i then one more application of (14.91) gives 

-7 7 pks(/0^) < 7pKs\P°] + ZC3T 2 < £[p ] +327r+2||/3 ||z,i(iK=,log(e+|xP)da ! ) + ZC 3 t 2 := Tq + Fit 2 , 
where (|4.11l) was used. ■ 

We are now ready to prove the Holder continuity estimate. 

5.3 LEMMA (Holder continuity). There is a positive constant T 2 depending only on the initial data and the 
regularization mollifier 7 such that for each r < f * and each k £ N, such that for all t > s > 0, 

W 2 (p T {t),p T {s)) < T 2 (t~s) 1/2 . 

Proof: Let j be such that (j — l)r < s < jr and let I be such that It < t < (I + l)r. By the geodesic property of 
McCann's displacement interpolation, 

W 2 (p T (s),pi) = 3 -^^W 2 (p>-\pi) and W 2 (p T (t) 7 pi) = t —W 2 (p i T ,p e + 1 ) . 

T T 

By Lemma I5T21 and the Log-HLS inequality (|4.10|) . .FpksIPt] ~ -^pksIPt -1 ] < -^b + F\t 2 - 87r(log8 - 1), and thus, 
by plugging into (|5.10|) . we get 

WKp^p*- 1 ) < 2r [T +T1T 2 - 8^(log8 - 1)] + 2ZC 3 T 3 2~ k 

< t \2T + 2T X - 16vr(log8 - 1) + 2ZC 3 ] := jf r , (5.12) 
since r < 1 and k £ N. Therefore, we deduce 

W 2 (p T ( S ),pi)< ] -l^F 2 ^ and W 2 (p T (t),pi)<— -T^ . 

T T 

Adding these two estimates and using the concavity of square root, 



W 2 (p T ( S ), ft) + W 2 (p T (t), p e T ) < T 2y /(t -*)-(£- j)r . (5.13) 

Next, by the triangle inequality, the Cauchy- Schwartz inequality, and (|5.10p and proceeding as in (|5.12[) . we finally 
conclude 

\ 1/2 

W 2 (/ T ,tfr 1 ) < J2W^Pr' X ) < V- jf' 2 ( E ( 2T [^PKSI^ 1 ] --FpKS&fl] +2ZC 3 T 3 T 
k=j \k=j 

<[(£-j)r] 1/2 (2 [JpKsIpr 1 ] -JvksIpR] +2ZC 3 t 2 ) 1/2 < T 2 {(I-j)t} 1 / 2 . 
Adding this to the estimate in (|5.13|) , and using the subadditivity of the square root concludes the proof. 
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5.3 Weak compactness 

In this subsection, we will show the compactness of the sequence of interpolating curves. We cannot proceed as 
usually done, for instance in [2], since we want to show that the limiting curves are not only measures but rather 
densities for each time and also since our densities have infinite second moment. The idea is to show a compactness 
in metrics W p , 1 < p < 2, less strong than W% and then, pass to the limit p /~ 2 the Holder continuity property. 
The next lemma shows a compactness property of the sets {p T (t) ■ < r < f*} for each fixed t. 

5.4 LEMMA (Uniform integrability at fixed t). There is a finite and computable constant F3 depending only on 
po and for any fixed 1 < p < 2 so that for all r < f * 

£[p T {t)\ = I p T (t, x) log/9 T (i, x) dx < F3 and I \x\ p p T (t, x) dx < F3 ■ 
Jr 2 ' Jr 2 

Proof: By the uniform control that we have on T-L\, the moment bound is immediate from Lemma 11.111 By the 
uniform control of J-pks and H\ in (|4.26p . Lemma [5.21 and by the first concentration control Theorem 11.91 we 
conclude that 

71 / p*log + p*da;< F + F lT 2 + C CC f (5.14) 
Jr 2 

where < 71 < 1, uniformly in k. Finally, using the bound of the absolute first moment together with (|5.14[) . we 
conclude that £[p^} is bounded uniformly in k. Then, by the displacement convexity of £ , this bound extends to 
p T (t) for all t > 0, as explained at the beginning of this subsection. ■ 



It follows immediately from (|2.5|) that Lemma 15.31 remains true if W2 there is replaced by any of the weaker 
metrics W p , 1 < p < 2. The following characterization of the convergence in W p metrics in [HJ Chapter 9]: 
convergence of the absolute p-moment plus the weak-* convergence as measures of a sequence of densities {p n }neN 
towards p is equivalent to W p convergence; implies the following compactness result. 

5.5 LEMMA (Compactness for the W p metric). For any M > 0, let K be a subset of the set of densities p of 
mass M that is uniformly integrable, and such that {\x\ p p{x) : p £ K} is also uniformly integrable. Suppose also 
that K is closed in the L x -weak topology. Then K is compact in the W p metric. 

Proof: Let {pn} n eN be any sequence in K. Since K is uniformly integrable and weakly closed in L 1 , the Dunford- 
Pettis Theorem provides us with a p € K and a subsequence {/?n fc }fceN such that limfc_ > . 00 p nk — p weakly in L 1 
and thus, weakly-* as measures. It is trivial to check that weak-L 1 convergence plus the uniform integrability of 
{|x| p p(a;) : p € K} implies that 

lim / \x\ p p nk (x) dx — / \x\ p p(x) dx . 
k ->°° Jr 2 Jr 2 

The characterization of W p convergence mentioned above then implies limfe_ i . 00 W p (p n( . , p) = 0. I 

5.6 THEOREM (Convergence asT-> 0). Given T > and any 1 < p < 2, define (M^,W P ) to be the metric 
space in which Ai^ is the set of densities on R 2 satisfying £[p] < and 

[ \x\ p+5 p{x)dx<F 3 (l+T^ p+s '>/ 2 ), 
Jr 2 

with p < p + S < 2. Then there is a function p on [0,oo) with values in the set of densities of mass 8tt such that 
for allT > 0, the restriction of p to [0,T] is continuous in (A^^,W P ) 7 and there is a sequence {r( n ^}„ e N such that 
for all T > 0, 

lim 

n— >oc 

Moreover the sequence {T'™'}„ g N can be chosen independently of p, i.e., such that the convergence property (|5.15j) 
holds for all 1 < p < 2. Furthermore, as a consequence for all t > s > and all 1 < p < 2: 

W p (p(t),p(s)) <F 2 (t~s) 1/2 . (5.16) 



max W„(p T 

0<t<T F 



lim 



max W (p_(,, 

0<t<T ' 



(t), P (t)) 



(5.15) 
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Proof: For each T > 0, {M. T , W p ) is a compact metric space as a consequence of Lemma [5.51 By Lemma | 
for each t <T, and each < r < f*, the restriction of p T to [0,T] takes values in (A^, W p ). Next, by (]2 ,5f) and 
Lemma [531 the set of these functions for < r < r* is uniformly equicontinuous into (A'ffi, W p ). Thus by the the 
Arzela-Ascoli Theorem, we can select a uniformly convergent sequence. 

Now a simple diagonal sequence argument concludes the construction of p and proof of (|5.15[) . Apply the above 
for T — 1 to get the initial sequence. Now take T = 2 and choose a subsequence of the first sequence, and so forth. 
For the piecewise-interpolation sequence, note that W2(/5 T (t), p T (t)) <Cyfr by standard properties of displacement 
interpolation and (|5.12[) in Lemma 15.31 Therefore, the limits t^ 1 ' — > of both time interpolations are the same. 
Note that this argument does not invoke any uniqueness of the limit. 

Another simple diagonal sequence argument shows that the sequence can be made independent of p. Take a 
sequence of increasing exponents {p n }nen /* 2 and {5 n } nS N \ with 1 < p n < p n + 6 n < 2 and apply the diagonal 
sequence argument to the constructed sequences for each p n . Also, take into account that the sequence of distances 
W p (p T („) (t),p(t)) is increasing in p. 

The last part of the claim (|5.16l) follows directly from (I5.15|) and Lemma 15.31 for all 1 < p < 2. Since the 
constant Ti obtained in Lemma T5. 31 does not depend on 1 < p < 2, then we conclude (j5.16[) for p = 2. ■ 

5.4 IP regularity 

Our goal in this section is to prove: 

5.7 THEOREM (X p -regularity) . For each finite a > and p > 1, there is a finite and computable constant C p 
depending only on a, p and po such that whenever r < a, 

p^(t, x) dx < C p for all t > a . 

M 2 

Proof: From (l4~27l) . we deduce that for each m, H\[p™} + V[p™\ < Hxip™- 1 } + ^2~"V 2 , proceeding in the 
same way that we did in deducing (|5.6j) . except this time we do not discard the dissipation term. Let n > k be 
positive integers. Since Q m < C Po for all m, summing from m = n — k to n yields 

Hx[pW]+ J2 rV[p'?]<<Hx[p n T - k - 1 ] + ^ . (5.17) 

m=n—k 

Then since < rl\[p^) < C Po for all k, using r < 1 and dividing by kr, we may simplify this to 

I E W]<^. (5-18) 

m—n—k 

We now choose k to be the greatest integer less than or equal to a/r, and of course suppose that n > k. Since 
kr < a < (k + l)r, and k > 1, a/2 < fcr, and then the fact that averages dominate minima yields the conclusion 
that for some positive integer m with r < a, 

4C 
a 

Then since H\[p™] < C Po and FpksIpT] < +°° we have from Theorem 11.101 that 



/ ivcc) 1 / 4 ! 2 dx < 



jftCpo CCCD 

072 72 



Recall the GNS inequality 

11/112 <^iiv/nr 4 n/ii4 
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valid for locally integrable functions / in R 2 and q > 4. Applying this with q = 4p and / = (p™) 1 / 4 , we obtain 



\\p?F p < 



B, 



47rC 0n a 



po , '-'CCD 



&72 



+ 



72 



2(p-l) 



8tt 



1/4 



(5.19) 



Thus we have an a-priori bound on for some m with n — k < m < n. We now apply Lemma |4.16[ in which 

the constant A\ and A 2 in (|5.20|) are defined, to conclude that 



/ ( P i n) - Kf + dx < f (p™ - K)\ dx + A 1 kT + A 2 Y / rV\p> T \ . 
Jr 2 Jr 2 ~L 



(5.20) 



We bound the first term on the right in (|5.20l) using (|5.19[l . the second using the fact that fcr < a, and the third 
using ([5J8|) . The result, 



(pW - K) p + dx < 



B 



ip 



\ a-f 2 72 / 



2(p-l) 



8tt 



1/4 



A ia + 2A 2 C, 



Pu 



uniformly for all n such that nr > a. Note that the bound depends only on a, p. Since ||pi ||p < || ipr ~ 

have the same type of bound on ||pi ||p> uniformly for all n such that nr > a. By 

the displacement convexity of 

pi-^ / p v {x) dx 



for p > 1, this bound immediately extends to p T {t) for all t > 



5.5 Verification that p = lim^o p T is a solution of the PKS system 

Let rW, p T<n) and p be given as in Theorem 15.61 Our goal in this subsection is to prove that p is a weak solution 
of the PKS system as given in item (1.5.2) of the Definition 11.51 

5.8 LEMMA (p is a weak solution of the PKS system). Let t^ 1 \ p r <n) and p be given as in Theorem lbM Then 
for all smooth and compactly supported test functions £ and all t 2 > ti > 0, 



/ 

Jr- 



C(x) [p(t 2 ,x) - p{ti,x)} dx=- — 

4tt j t 



, , , x (g - y) ■ (vOf) - vc(y)) , , 

p(s, a;) p(s, y) — — dy dx 



\x - y\ 



+ I I AC( X ) p( s i x ) da; ds , 



Proof: In order to prove Lemma 15.81 we first remind the reader an analog for the functions p T : 

5.9 LEMMA (Approximate weak solutions of the PKS system). For < r < 1, define e[t,r) — for t 6 
((fc — l)r, kr], and all integers k > 1. Then for all smooth and compactly supported test functions £ and all 
t 2 >h> 0, 



C(x) [p T (t 2 ,x) - p T (t 1 ,x)]dx 



1 



2 xB 2 



p T (s,x) p T (s,y)VG e ( t . T) (x - y) ■ (VC(x) - VC(y)) dydx 



+ 



AC(x) p r (s,a;) da; ds + 0(r 1/2 ) . 



The proof of the previous Lemma follows the same lines as in [5] Theorem 3.4] and we skip it here for the sake 
of conciseness. The interested reader can see its proof in the preprint version of this paper. 
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As r — > along the sequence {r(™)} ne N, lim„^.oo Wi[p T (*j (t), p(t)] = uniformly on [0, T] for any finite T. 
Hence by the Kantorovich-Rubinstein Theorem, 



lim / A£(x) p r (n) it, x) dx = / A£(x) x) dx , 

™^°° JR 2 Jr 2 

uniformly on [0,T]. The interaction term can be rewritten as 

/ [VG e(t)T („)) * p rM [t)] (x) ■ V((x) p rM (t, x) dx 
Jr 2 

= / [VG * 7e(t.r(")) * P r M (*)] (x) • [7 e (t,r<")) * (Pt<»> (*) VQ] (x) dx 

1/7 ( w, w w, v (x - y) ■ (VC(x) - VC(y)) , , 

= IT Ue*PTM){t,x){-j e *p T (n)){t,y) : 12 dxdy 

47tJJ K 2 xR 2 |x-y| z 

+ / [VG * 7 e(tiT (n )) * p T in) (t)] (x) • [7 e (t,r<">) * (PrC») (*) VC) - (7 £ (t.r(">) * p r M (*))VC] (x) dx 

JR 2 

:=h+h- (5.21) 

As {p r <n)(t)}neN converges weakly in L 1 (K 2 ) towards p(t) as n — > oo, so does {7e(t, T 0»)) * /0 r (») (t)}neN- We then 
deduce that {(7 e (t >T («)) * p T (n)(t)) (g> (7 e (t )T c»)) * P r («) converges weakly in L 1 (R 2 x R 2 ) towards ® p(t) 
when n — > oo, see [5] Lemma 2.3]. As a consequence we can pass to the limit in the first term in the right-hand-side 
of (|E2T|) to obtain 



hm h = -— / / x) y) : (2 dx dy 

47rjj R 2 XR 2 \x-y\ 



We must now show that la disappears in the limit. We can estimate /2 using 

|7e(t,T<«>) * (PrW (*) VC)-(7e(t,r(«)) * Pr« (*))VC|(x) < / 7e(t,r<»))(» " 2/) l V C(Z/) - V((x)| p r („, (t, y) dy 

JR 2 

< Cq / 7e(t,r(»))( a; - V)\ x ~ V\ ArC»>(*>2/) d 2/ = ((7e(t,T<»>) |x|) * p T M (t))(x) . 
JR 2 

By the HLS inequality, (|H)j) . 

|/2| < ^^||7 e (t,r<")) * Pr<") WIU/3 Il7e(t,x(«)) * (PrW (*) V ~ (7 e (t,r(">) * /M") (*)) VCIU/3- 

Then by similar arguments similar to those used to prove Lemma |4.13[ we get 



\h I < 4 Ghls G c II p rM (t) || 4/3 II 7 |x| || 4/3 yfittyrtj ■ 



In case t\ > 0, estimating ||/9 T (»») (0 II 4/3 using Theorem 15. 7\ we obtain the result. If t% — we can use instead (|4~29|) 
and Holder's inequality to obtain 



\h\ < 4G HLS G c (Stt)^ V^hkllU/s^^ = 0((r M ) 7 / 6 ) , 
where (|5.4I) was used. ■ 

5.6 Strong Compactness 

At this point we have shown that the limit p = lim n _>oo p T (n) posesses the properties (1.5.1) and (1.5.2) in Defini- 
tion Q~3] of properly dissipative weak solutions. In this subsection, we show that (1.5.3) is also satisfied. This will 



BCC August 1, 2011 



45 



complete the proof of the existence of properly dissipative solutions in Theorem 11.61 Choosing n — k in (|5.17[) we 
obtain, for all k € N and all r = r( n) that 

Hx\ffi + £ rV[p™} < n x [ Po ] + ^fr 2 . 

m— 1 

Thus, using the discrete time interpolation p T , we have that for any T > and the positive integer N such that 
Nt <T < (N + 1)T, 

Ux[?t{T)] + / V[p T {t)] dt < Hx[p } + -^r 2 . (5.22) 
Jo 4 

We emphasize that the use of the piecewise constant interpolation is essential at this point since the functional 

V[p] is not displacement convex. 

Note that the LP bounds deduced in Theorem 15 . 71 apply to p r as well as to p T . To make full use of these bounds, 

we choose any fixed a > 0, and then for all r < a, we weaken the bound in (|5.22l) by increasing the lower limit of 

integration in t to a. Also writing b :— Nt, this yields 

Hx[pr(T)}+8 [ [ \Mp l J\t lX )\ 2 dt- f ( pl/ 2 (t,x)dxdt<n x [po} + ^r 2 . 

Ja JM 2 Ja JK, 2 4 

It is legitimate to express T>[p T ] as the difference of two integrals since Theorem 15.71 tells us the f?J 2 is integrable 
over [a, T] x R 2 . We now show that passing to a further subsequence of {t(")}„ g n, we may arrange that for all 
< a < b < oo, along this subsequence, 

lim / / \p T(n) (t,x) - p(t,x)\ 3/2 dxdt = (5.23) 



and 

lim p T(n) (t,x) = p(t,x) 

n— >oo 

for almost every (t, x) £ [a, b] xl 2 . The following strong compactness theorem, leading to the existence of almost 
everywhere convergent subsequences, is the key: 

5.10 THEOREM (Strong compactness for p 2 (t,x)). Let < a < b < oo be given. The family of functions 

{ Ma.b]{t)Pr(t,x) : T < T* } 

is precompact in the strong topology in L 2 (R 2 x (a, b)). 

Proof: By the Kolmogorov Compactness Theorem [23l [28] , also known as the Frechet-Kolmogorov theorem j£l 
Corollary IV. 26], a set Q of functions g £ L 2 ((a, b) x R 2 ) is strongly precompact if and only if: 

(Kl) Q is uniformly bounded in L 2 ((a,b) x M 2 ). 

(K2) For each e > 0, there is an R > so that for all g G Q, 



[ [ \g(t,x)\ 2 dxdt < e. 

Ja J\x\>R 

(K3) For each e > 0, there is a S > so that for all g € G, whenever y £R 2 satisfy \y\ < 5, 

\g{t,x~y)-g(t,x)\ 2 dxdt < e , (5.24) 



! x(a,b) 

and whenever c G R, and [a, b] C (a, b) with (a — c,b — c) C (a, b) 



\g(t- c,x) - g(t,x)\ 2 dxdt < e . (5.25) 

! x(a,b) 
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Step 1.- Uniform integrability for p 2 (t, x): First, note that for each r, p 2 = p¥ 2 pr ■ Therefore, by the Cauchy- 
Schwartz inequality 

\x\ 1/2 p 2 T (t,x) dx < ( [ \x\p T (t,x)dx) ([ p 3 T (t,x)dx > 

\JR 2 J \JR 2 J 

Recall the first moment is controlled by T-L\ in Lemma and that H\[p T (t)} < C Po . Also, Theorem 15.71 give us 
a bound on ||/5 T (i)||3 uniformly in t £ [a,b] for all sufficiently small r. Thus, there is a constant C depending only 
on a and b so that for all sufficiently small t, 

b 

\x\ 1/2 p~ 2 T (t,x) dxdt < C 



Even more simply, by Theorem 15.71 we have a constant C depending only on a and b so that for all sufficiently 
small t, such that 

/ / p 2 T (t,x) dxdt < C and / / p%{t, x) dx dt < C . 

Ja JR 2 J a JR 2 

The uniform integrability is an immediate consequence of these estimates, giving the first two conditions (Kl) and 
(K2) of the Kolmogorov compactness theorem. 

Step 2.- Spatial translations (|5.24[) : Writing p T (t,x) = f (t,x), we have 



Pr(t, x-y)- p T (t, x) = -A\y\ ■ I f (t,x- sy)Vf(t, x ~ sb) ds , 

J o 

and then by Minkowski's inequality, and then Holder's inequality 

\p T (t,x-y)-p T (t,x)\ 2 dx < 4|y| ^J 6 (t,x)dx^ ' f(t,x)\ 2 dx^j ' 

= 4|y| ^0 2 (t,x)d x y 2 |VpV4|2 d ,y /2 . 
By our uniform L p bounds, there is a constant C independent of r so that 

1 /2 

f \p T (t,x-y)-p T {t,x)\ 2 dx< \y\c( f |Vpy 4 | 2 dx) , 

JR 2 \JR 2 / 

and hence 

b f b \ 1 I 2 

f f \p T (t,x-y) - p T (t,x)\ 2 dxdt < C\y\b 1/2 ( [ f \V p 1 ^ 2 dx dt) . 

Ja JR 2 \Ja JR 2 ) 

Since the integral on the right hand side is bounded uniformly in r, this gives us ()5.24|) . 

Step 3.- Temporal translations (|5.25l) : We do not have any estimates on time derivatives, so we cannot obtain 
the bound on temporal translations in such a simple manner as we have for the spatial translations. What we do 
have from Lemma 15.31 is a finite constant C so that 

W a (p(t-c,.),p(t,-))<CyS (5.26) 

holds uniformly in r < f * and in t > c. 

We now use an interpolation argument based on an idea of Otto, see |36l Subsection 3.5], to combine this with 
the spatial regularity provided by the square integrability of V p¥ 4 (t,x) in R 2 x R. Our task would be very much 
simpler if we had, for almost every t, a uniform bound on L 2 | V/5^ 4 (t, x)\ 2 dx, but don't for any fixed t. 

The interpolation bound we need is provided by the following theorem whose proof will be provided after we 
conclude the proof of Theorem l5.10l 
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5.11 THEOREM (Interpolation bound). Let cto and a\ be two densities of mass M in R 2 such that for some 
p > 2, \\o~o\\p^\ , ||ci||p+i < K, Suppose also that o\J^ and o~\^ have square integrable distributional gradients. 
Then 

ho-vi\\l < (||V( CT0 ) 1 / 4 || 2 + ||V( CT1 ) 1/4 ||2) (2^+2^K)(W 2 (a ,a 1 ))^/^ 

We now apply Theorem 15 . 1 1 1 with p T (t, •) in place of Co and p T (t — c, •) in place of o\. We have a uniform bound 
on K in this case, for any p < oo, and we also have the bound (I5.26|) . so we obtain a finite constant C so that 



|c| & p+ 4 + |c| 4 ?+ 2 



lft.(t,0-M*-c,-)rd»< 

o ( (| r2 vp-y 4 (t, -)i 2 da) 1/2 + (/ r2 v P -y 4 (* - c oi 2 dx) 1/2 

Now integrating both sides over [5, b], and using the Cauchy-Schwartz inequality, we obtain 

/ \Pr(t,-) -p T (t-C,-)\ 2 dxdt < C(AT 1 / 2 \c\^ +T\ C \^) , 

JR 2 x[d,b] \ ' 

where 

V2 / „ T „ \ 1/2 



A := 



f / T / V#/ 4 (t, -)| 2 dx + ( f f Vp¥*(t - c, -)\ 2 dx dt) 

\Jto JR 2 J \Jto JR 2 ) 



Our results so far give us a bound on A that is uniform in r < f* , and thus, choosing p — 2, the proof of (|5.25[) is 
complete. ■ 

It remains to prove Theorem l5.11l Before beginning the proof itself, we explain the argument in [551 Subsection 
3.5] that is the basis of the proof. Let ctq and u\ be two uniformly bounded densities of mass M in R 2 . Let dLT 
be the optimal coupling of o~o and o~\. For < s < 1, define a s to be the displacement interpolant between ao and 
el . That is, for any bounded continuous function tp in R 2 , 



ip(z)a s (z) dz = I <p((l - s)x + sy) dU(x,y) . 

JR 2 xR 2 

As is well known (see for example |37[ Lemma 3]), we have 

Ikslloo < max{ \\a \\oo , Ikilloo } • (5-27) 

Suppose furthermore that both Co and o~\ have a square integrable gradient. Then we have Otto's interpolation 
estimate: 

\\o~Q — 0"1 1| 2 < (max{ ||(7o||oo > ||ci||oo }) 1/2 (||Vao|| 2 + ||Vffi|| 2 )W2 (cto, ci) . (5.28) 
To see this, note that 

Iko-crill 2 = / ([o-a(x)-ai(x)}-[a Q (y)-a 1 {y)})dU(x,y) 

[Voi(x + s(y — x)) — Vcro(a; + s(y — x))} ■ (x — y) dIT(x, y) ds . 

;R 2 

Now apply the Cauchy-Schwartz inequality. By (|5.27|) . we have that for each s, 

/ |Va (a; + s(t/-a;))| 2 dn(x,2/) = / |VCT (a;)| 2 cr s (a;) dz < (max{ ||o- ||oo , IMIoc }) ||Va ||^ , 

JR 2 xR 2 JR 2 
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and likewise for <j\. Since J R2xR 2 I* — y\ 2 dll(x, y) = W|(cro, we obtain (|5.28l) . 

We would like to apply this sort of argument with p T (t, •) in place of Co and p T (t — c, •) in place of o~\. We 
cannot do this directly since p T is not bounded. However, by the results of the previous section, for to > 0, we have 
an upper bound on ||/5 r (i, -)||p for all finite p that is uniform in t > t . The next lemma allows us to approximate 
p T (t, •) and p T (t — c, •) by uniformly bounded densities without significantly increasing the 2-Wasserstein distance 
between them. 

5.12 LEMMA (Approximation by bounded densities). Let a he any density of mass M on R 2 . Suppose that for 
some p > 1, a belongs to L P+1 (R 2 ) with ||o'||p+ 1 < K . Then for all A > 0, there exists density of mass M, a such 
that Halloo < 2A, 

W 2 (o» < 2K\- p - 1/2 , (5.29) 
||cr - oF|| 2 < 2K 1 ^\- ( - p - 1 ^ 2 , (5.30) 

and finally, ||5|| 2 < 2 1 1 cr 1 1 2 . 

Proof: Let E\ = { x : a(x) > A }. Then by Chebychev's inequality, 

/ a{x)dx< i / aP+ 1 (x)dx<^- . 

Now for any h > and any integers m and n, define C mj „ to be the square 

C m ,n — {x = {x\, X2) G R 2 : mh < x\ < (m + l)h , nh < X2 < (n + l)ft } . 

Define M m ^ n by 



(5.31) 



Next define 



0~ m ,n '■— l-E x nC mn V 



a(x) dx 



E x r\C„ 



and 



There are both densities of mass M m ^ n supported in C mi „. Let T mt „ be any map from C m>n to C TOi „ such that 
Tm,n4f z(T m.n = ^m.n- For instance, one may use the one associated to the optimal coupling. Define T : R 2 — > R 2 by 



T(x) 



^ lc m , n ( x ) T ?n,n(x) X E Ex 

xiE x 



Notice that only one term in the sum is non-zero. By construction, 

r#a = (i-i B >+ ^tF 1 ^ 



By (|QT1) . M m , n < if A _p . We now specify ft := \-(p+i)/2 _ Then M m ^ n h~ 2 < A, and hence with this choice 
of ft, T#ct is uniformly bounded above by 2A. We now define a := T#cr. It remains to verify that a has all of the 
properties claimed in the lemma. 

First of all, notice that \T(x) — x\ 2 = on the complement of E\, and on E\ is is bounded by 2ft 2 , since each 
T m ,n maps C m>n into itself. Therefore, 

W2(ff, 5=) < / l^(x) - x\ 2 a{x) dz < 2ft 2 f tr(x) dx < 2K\~ ( - P+1 ^K\~ P . 

JR. 2 J E x 

This proves (|5.29p . Next, by (|5.31[) and Jensen's inequality, 



~ 1 \ " Mva,n -. 

0-0= ±E x (r- 2^ —f^-^m.n 



with 



h 2 



< UeM 



(5.32) 
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Furthermore, by Holder's inequality with the dual indices p and p' , 

/ a 2 l Ex dx = [ a 2 - 1 !'?' o 1 lv'\ Ex &x < ([ a p+1 dx] ^ ( [ a dx) ^ 
Jr 2 Jr 2 \Jr 2 J \Je x J 

< (K) 1/p (K\- p ) 1/p ' < K\ l ~ p . (5.33) 

This proves (fOOj) . For the final part, note that by (|Q2]) . ||ct|| 2 < ||(1 — l jEA )cr||2 + ||l£; A cr||2 < 2||ct|| 2 . This completes 
the proof of the lemma. I 

Proof of Theorem 15. lit For any A > 0, let ao and a\ be the bounded approximation of cto and <7\ provided by 
Lemma T5. 121 Also, define 

up (a:) := min{cri(a;) , A } and a\[x) := min{<7i(x) , A } . 

Starting from the identity ||<7o — ai\\ 2 = {cq — a\,oo)h 2 — (&o — cij om,)i, 2 ) and adding and subtracting repeatedly, 
and using the Cauchy- Schwartz inequality, we obtain 

Iko — 01II2 ^ (go ~ ^i.^o)l 2 - (go - oj_,&i)l 2 

+ (Ikolla + 1 1 o"i ||2)(|| - co 1 1 2 + 1 1 o-i - CT1II2) 

+ (Hffolla + llffilla) [(Iko-^olb + lki-^ilb)] • (5.34) 

The heart of the matter is the estimation of (oq — aj_, <xo)l 2 — (go — &i)l 2 - Let LT be the optimal coupling of 
<7o and (7i, and for < s < 1, let a s be displacement interpolant between ao and <7i. By what has been explained 
above, a s is uniformly bounded above by 2A. Then, by the Cauchy-Schwarz inequality once more, 

(ao- £1, (j )l 2 - oj_,vi)l 2 = / (\go -£i|(x) - [£p_ -£i](y)) dll(a;,t/) 

JR 2 xR 2 

= 11 ([Vaa-Va{\((l-s)y + sx)-(x-y)dn(x,y) 

JQ JR 2 xR 2 

< (2A) 1 / 2 (||Voo||2 + ||Wi|| 2 ) W 2 (5 : o,ffi) ■ 
Next, for j = 0, 1, let f 3 := (ctj) 1/4 - Then since Va 3 (x) = unless fj(x) < A 1 / 4 , ||Vct^|| 2 < 4A 3 / 4 || V/j|| 2 . Thus, 

(ao - 01,0-0) l 2 ~ (oo - cti,?i)l 2 < 2 5 / 2 A 5 / 4 (|| V(a ) 1/4 || 2 + || V^) 1 / 4 ^) W 2 (a ,5 1 ) . (5.35) 
By the triangle inequality and Lemma I5.12| 

W 2 (ff ,CTi) < W 2 (ct ,CT ) + W 2 (cr ,CTl) + W 2 ((7i,CTi) 

< 2KX-P- 1 / 2 + W 2 (a ,a 1 ) + 2K\- p - 1 ' 2 . 

Combining this with f)5 -35[) . we obtain 

(ao-ai,a ) L 2 - (ao -£i,2i) L > < 2 5 / 2 A 5 / 4 (||V(a ) 1/4 || 2 + HVfa) 1 / 4 ^) W 2 (a , f r 1 ) 

+ 2 5 / 2 A 5 / 4 (\\V(a )^% + WVia^/%) AKX^ 2 . 

At this point we specify 

A:=(W 2 (a , ( 7 1 ))- 2 /( 2 f+ 1 ) . (5.36) 

Wit this choice, we have 

(og-£i,3F ) l2 -(ao-ai,?i) l2 < (\Ma ) 1/4 h + WnVh) (2 5/2 + 2 9 / 2 ^)(W 2 (a , CTi ))(4p-3)/(4p+2) . 
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Finally, it remains to bound the terms in the last two lines of (|5.34p . Note that the same argument using 
Chebychev's inequality in (|5.31l) and then Holder's inequality in (|5.33l) yields \\<7j — cr^- 1 1 2 < K 1 / 2 \( 1 ~ p ' ) / 2 for 
j = 0, 1. Note also that by Holder's inequality once more, \\o-j\\ 2 < M( p ~ 1 ^ p K 1 l p . Therefore, the sum of the terms 
in the last two lines of (j5~34| is bounded above by lQM^- 1 ^ p K ( ^ p+2 ^ 2p \-'- p - 1 ^ 2 . With the value of A specified in 
(|5.36l) . the contribution of these terms is 

Combining results, the proof is complete. I 



5.13 LEMMA (Further subsequence). There is a subsequence of the sequence {r'"'}„ e n, denoted with the same 
index, such that along this subsequence, (|5.23p is valid for each < a < b < oo and 

lim j5 t (t.) (t, x) = p(t, x) almost everywhere on (0, T) x R 2 

n— >oo 

for any T > 0. 

Proof: Let us consider any integer N > N a with N a T > 1. Applying Theorem 15.101 we get a subsequence of 
{r^jneN along which {p T ( n )} n eN converges to p strongly in L 2 ([l/N , T]xR 2 ) and almost surely on [1/N , T]xE 2 . 
Next, for N = N + l, choose a subsequence of {r^ n '} of the previous subsequence along which {p T ( n )} n eN converges 
to p strongly in L 2 ([1/(A D + 1),T] x R 2 ) and almost surely on [l/(N + 1),T] x R 2 , and so forth. We finish by an 
obvious diagonal sequence argument. ■ 



5.7 Entropy dissipation 

With the strong convergence results obtained in the previous subsection, we may now establish the entropy-entropy 
dissipation inequality: 

5.14 THEOREM (Entropy-entropy dissipation). For each T > the weak solution p of the PKS system that we 
have constructed for the initial data po satisfies 



J~pks[p{T)] < +oo and H x [p(T)} + [ V[p{t)] At < H x [p } 

Jo 



(5.37) 



Proof: First remind that J-pks is lower semi-continuous function for the weak-L 1 convergence, see [51 Lemma 3.1]. 
So that the bound in (15.37)) is a direct consequence of (|5.23p . Also by (|5.23p . it suffices to show that 

Hx\p(T)]<]imMHx\? r w(T)] , 



and that for any 0<a<6<T<oo 



V (p 



,V4 



dx dt < lim inf 



a Jm 2 



v (p% 



dx dt . 



(5.38) 



for a suitable sequence {T^} ne N, since the rest easily follows by a monotone convergence argument for taking a 
to and b to T. 

The first of these follows from the fact that T-L\[p] is a lower semi-continuous function on L 1 with respect to the 
Wi metric just by using the expression of H\[p] in (jl.lip . 

To see the second, we again make use of the almost everywhere convergence proved in the previous subsection. 
Let /„ := p 1 '^ and / := p 1 / 4 . Then the sequence of functions {/ n } n eN — > / in L 4 n L 6 ((a,T) x M 2 ) from 
Lemma T5.13I From (|5.22p . we have that the sequence {V/njngN is bounded in L 2 ((a,T) x R 2 ), therefore it has a 
weakly convergent subsequence denoted with the same index such that {V/„}„ e N — 1 cr weakly in L 2 ((a,T) x R 2 ). 
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Due to the strong convergence of the sequence {/ n } n6 n — > f in i 4 n L e ((a,T) x R 2 ), it is simple to identify the 
weak limit as a = V/. By standard properties of L 2 -weak convergence, we deduce that 



/ / |V/| 2 da; dt < liminf / / |V/ n | 2 da; dt 



which shows f|5 . 38|) . 



Proof of Theorems 1 1.6t As noted above, Theorem 15 . 141 provides the final step in the construction of the properly 
dissipative weak solutions. Theorems 15 . 71 and 15 . 141 provide the additional regularity properties (1.6.1) and (1.6.2). 
It remains to prove (1.6.3), the dissipation of J-pks- 

We now show that J-pKs[p{t)] < Fpks[p{s)] for all < s < t. Take r to be any element of the sequence {r(™)} ng N 
whose corresponding approximated solutions {p T ( n )} n £N converges to the constructed properly dissipative weak 
solution p. Let j be such that (j — l)r < s < jr and let I be such that It < t < (£+ l)r. Using (|5.11l) in Lemma I5T21 
and (14.91). wc deduce 



T PKS [Pr(t)] < T^ KS [pU < J~pKg [p 1 - 1 ] + ZC 3 r 2 2-t < [pi- 1 ] + ZC 3 r 2 £ 2" fe < ^ [p^ 1 ] + ZC 3 r 2 . 

k=j 

Using (|5.9j) and the Lemma 15.71 we can control the error term in the right-hand side by 

^p KS [Pr(t)] < T PKS [p T (s)] + Z\\p t (s)\\1 / / It 3 2-i +1 + ZC 3 t 2 < T PKS [p T (s)] + 2ZC 3 t 2 . 

Finally, the a-priori bounds uniform in t due to Lemmas 15.61 and 15.71 together with Lemma 15.161 allow us to pass 
to the limit r — > leading to our claim. I 

5.8 Large time asymptotics 

We start by identifying the large time asymptotics of the solutions in a time average sense. 

5.15 LEMMA (Time-averaged strong convergence). Let p be the properly dissipative weak solution of the PKS 
system that we have constructed. Then 

T+l 
T 



lira ( / / \p(t,x)-Qx(x)\dxdt ) =0 . (5.39) 



Proof: This follows by a standard entropy dissipation argument. Let {i n }„ e N +°° be an increasing diverging se- 
quence of times and consider a n (t, x) = p(t+t n , x), for < t < 1. By using the entropy dissipation inequality (|5.37[) 
which is true for all T > 0, we deduce that 



/>oo /> 1 

/ V[p{t)] dt < oo , and thus, lim / T>[a n (t)] dt = . 
Jo ™^°°Jo 



(5.40) 



Now, again using the uniform in time bounds (|5.37l) for the solution, J r pKs[c r n(0] < +oo, and H\[cr n (t)} < 
H\[po] < C Po , the concentration control inequality in Lemma ll. Ill and the Gagliardo-Nirenberg-Sobolev inequality 
in Lemma [L2l we deduce 



i ft a^d X dt< H V[a^] 2 dxdt<^ [ 1 V[a n ]dt+^<^C po+ ^. (V41, 
o Jo in 2 Jo Jr 2 72 Jo 



72 72 72 



Moreover, due to Theorem 15 .71 we deduce 



sup sup / o*(t,x) dx < C p , (5-42) 

nGN 0<t<l , 
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for all 1 < p < oo. Note that the sequence {a n }nefi satisfies the equicontinuity property f|5 . 1 6[) in Theorem[ 

Summarizing, the sequence {<7„} n6 N has the same properties f|5-41[) and (|5.42l) as the sequence of approximate 
solutions we used in previous sections to construct the solution in Theorem 11.61 Proceeding as in Subsections 
5.3-5.6, we deduce the existence of a subsequence, denoted with the same index, such that {<T n } n gN converges 
towards p^ with the same convergence properties as in previous Subsections 5.3-5.6. Here, p^ is a weak solution 
of (jl.ip on the time interval (0, 1) in the sense of (1.5.2) in Definition [L5j In particular, {<r„}„ e N converges to p^ 

1/2 i 

in the metric space (A4/ , Wi), with the notation of Thcorem l5.61 giving 



lim 

n— >oo 



max Wi^^poo^)) 

0<t<l 







Moreover, repeating the arguments in Theorems 15. 141 and 15.131 we get 



sup Hxlpooit)] < C Po and liminf 



0<t<l 



o Jr 2 



VK) 



1/4 



d£ dt > 



o Jr 2 



Va 



1/4 



dx dt , 



(5.43) 



(5.44) 



and 



lim / / al /2 (t,x) dx dt = / / p 3 J, 2 (t, x) dx dt . 

Furthermore, Theorem 15.131 implies the almost everywhere convergence in (0, 1) x R 2 of {<7„} ne N towards poo, that 
together with f|5 -43[) implies that 



lim 



\a n (t,x) - poo(t,x)\ dxdt) = 



(5.45) 



Now, let us identify the limit function p^, passing to the limit using (|5.40p . we obtain 



Vp 



1/4 

oo 



V 2 ) dxdt = 0, 



which means that Poo(t) is a minimizer to the Gagliardo-Nirenberg-Sobolev inequality for all t £ (0,1), sec 
Lemma 11.21 and thus that there exists X(t) such that Poo(t) — Q\(t) where q\ is the family of the minimizers 
of the Gagliardo-Nirenberg-Sobolev inequality, see Lemma [1.21 Due to (|5.44[) then - Ha[£ , a(*)] < 00 1 we conclude 
that X(t) = A since H\[q^} = +oo for \i ^ A. Therefore, Poo{t) = q\ that together with (|5.45j) implies ()5.39|) . ■ 



We now will take advantage of the other Lyapunov functional, we shall prove that limt_>oo ■T'PKS [/>(£)] = 
^"pksIpa]- In doing this, we shall make essential use of the monotonicity of J 7 pks[p(^)]- Let us introduce for 
any C > and S > the set 



C,5 



\p £ L\{^ 2 ) : [ p{x) dx = 8tt , / \x\'p(x) dx < C and / p 1+S (x) dx < c\ 
I Js. 2 Jr. 2 Jr 2 J 



5.16 LEMMA (Convergence for .Fpks)- Given any sequence {p n }neN in Sc.s there is a p £ Sc,s and a subsequence 
{Pn k }ken su ch that 

lim llPnfc - Pill = and lim -Tries [PnjJ = FpKs[p] ■ 

Proof: Choose any < 6' < min(i5, 1) so that 25' / (1 — 5') < 5. By uniform integrability arguments such as we have 
made above, see Step 1 of Theorem 15.101 we can find a subsequence (denoted with the same index) along which 
{Pn}neN is weakly convergent in L 1+s (K 2 ) and along which {p]^ s }„gn is weakly convergent in L 1 . It follows as 
in Subsection 5.6 that {p n }n£N is strongly convergent in L 1+s (M 2 ), and passing to a further subsequence, we may 
suppose it is also almost everywhere convergent, and strongly convergent in ^(R 2 ). Let p denote the limit. By 
Fatou's Lemma, p £ Sc,s- 
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Since for t > 1, t\ogt < (t s - 1)/S', we have for p > 1, plogp < (l/8')p 1+6 ' and for p < 1, plog(l/p) < 
{l/S^p 1 - 5 ' . Since for e = 26', we get 



p 1 - 5 '(i + H 2 ) e d a; = / p^'ii + ixff'ii + ixfy'-dx 

Js. 2 

<(f pa + M 2 ) 2e/(i - 5,) v / (l+N^d^y , 

our choice of (5' gives the uniform integrability of {p n ~ s }neN- Then, by what we have said above, {p n log p n }ner 
is uniformly integrable, and hence 

lim / p n log p n dx = / plogp dx . 



The convergence of the positive part of the interaction potential is straightforward, due to the uniform bound of 
{NVnlneN in i 1 (M 2 ) and a dominated convergence argument. Concerning the negative part, it follows by Young's 
inequality for convolutions using the convergence of {p n }neN -> P in L 1+S '(R 2 ) and the fact that log_ |sc| € LP(R 2 ) 
for all 1 < p < oo. ■ 



5.17 LEMMA (Qualitative stability for .Fpks)- For any e > and C > 0, there exists 8(e,C) > so that if 
p G Sc,s, then 

•FpKsM < 87r(-l + log(87r)) + S(e,C) =>> e,Ji < e for some /i > , 

and /or any R > 0, t/iere exists 6(e, C, R) > suc/i i/ia£ 

1/2 

l(Vp-^)(^)| 2 d^] 

'{|*|<R} 



TpksW <87r(-l + log(87r))+ ( J(e,C ) fl) (7 KVp " V^)( x )\ 2 dx ^) ^ e for some /i > . 

W-fkl<R> / 



Proof: Given C,R > fixed, suppose not. Then for some e > 0, there is a sequence {p n }n&i in Sc,s such that 
lim^oo JpKs[Pn] = 8vr(-l + log(87r)) but 

inf \\p n - QpWi > e . 



However, by Lemma [5.161 there is a subsequence, still indexed by n, converging strongly in L 1 (K 2 ) to p € Sc. 8, 
such that 

8?r(-l + log(87r)) = lim 7pKs[Pn] = -^pksW • 

n— too 

By the cases of equality in the Log-HLS inequality, p = p M for some p. This is a contradiction. The second part is 
proved the same way, using the uniform integrability of the y/p, p € Sc,s on {\x\ < R}. I 

Proof of Theorems ll.7t Recall that J-pKs[p(t)} < ^pks[/o(s)] for all < s < t. We now apply this monotonicity 
to improve our large time asymptotic result. 

By (|5.39[) in Lemma T5. 151 there is a sequence of times {£ n }n6N /* co such that lim n _j.oo \\p(t n ) — Q\\\i — 0. By 
our regularity results in Lemmas 15.61 and 15. 7[ {p(t n )} ne ?j C Sc,s for some < C, 6 < oo. Then by Lemma [5.161 
there is a subsequence, still indexed by n, such that lim n _>oo J r pKs[p(^n)] = ^pksI^a]- By the monotonicity of 
J~PKs[p(t)] it follows that 

lim 7pks[p(*)] = ^pksIpa] = 8tt(-1 + log(87r)) . 

t—>oo 

Then by Lemma T5. 171 it follows that given R > there exists \x > such that for all sufficiently large t, 
||p(*)-e M ||i<e and llv^T" Vfllvi : = (7 KVp ~ vW*)! 2 da; l ^ e - 
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However, for any R > 0, by Minkowskii's inequality and (|5.37l) . 

1/2 



(f \V^-V^\ 2 e^ 1/2dx ) <([ Iv^T- Vp\ 2 0^ 1/2 ^] + 

\./{M<R} J \J{\x\<R} j 



<V8X(X + R 2 )\\^-^hM + ^C po . 

Since the left hand side diverges as R increases, uniformly for |p — A| > 5 > 0, we readily conclude that p, = A and 
]im t ^ oo \\p{t)-Qx\\i=0. ■ 



Appendix: proof of Lemma 14.131 

By Lemma 12.31 applied to uq — p and u\ = pa 



2xf' s (\xf) + L V " 



2 (p + <5) 3 / 2 



(V<p(x) - x) p dx - if«(p, p ) 



with /£(s) = 8A [8A + S(X + s) 2 ] 
it as 



-3/2 



Using (|4.20p . i.e. x — V<p(x) — r (Vc e — Vp/p) and expanding we can rewrite 



7/A.4H :• ?/A.//>«.j -rr ( |(I) + | (II) + 2 (III) + 2 (IV) ) - /O.p.p,,! 



where 



(I) 



iv P r 



dx, 



(II) := / d, 



(p + 5)3/2 > W y R2 ( p + 5)3/2 

(HI) := / f's(\x\ 2 )x- Vc e pdx and (IV) := - / /^(|x| 2 ) x • Vp dx , 

JR 2 JR 2 

We will keep the term (I) and we need to perform some integration-by-parts in the other terms: 
Control of (II): We can rewrite this term as 



(II) -2 / V 

JR 2 



p + 25 



Vc e dx . 



Integrating by parts on the ball of radius R and noticing that — Ac £ = 7 e * p * 7 e , we obtain 

p + 2<5 



\x\<R \VP + 3 



\x\<R. \VP + $ 
p + 25 

\x\=R Vp+~s 



(7 e * p * 7 e ) dx 

Vc e • n dcr 



< 2 / -y/p (7 e * p * 7 C ) dx + 4 \/# / 7e * p * 7 e dx 

J|a;|<_R J\x\<R 



p + 25 



IVcJ dtr 



< 



2/ ^(7^*7^x4-32^ 
+ 2 / (VP + 2V^) |Vce| da , 

J\x\=R 



where we used twice the estimate 



p + 25 
Vp~+7 



<Sp + 2V5 



(5.46) 
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Let us deal first with the second boundary term. By (|4.5[) in Lemma we have 

AVd f |Vc £ |da< W£(8C H Ls||7ll4/ 3 |IM7ll4/ 3 + 4+ ^ || 7 |U/, 

•J \x \ — R \ ~ ■ ' 1 



+ ^ Nil/all 14*) (5-47) 



>\x\=l 

:= 16ttVS{ J 7 + 1 + (%) , 

for all R > 0. To cope with the first boundary term, we observe that taking any natural N > 1, and considering 
( [ N [ VplVQldadr^ <n(N 2 -(N-l) 2 ) \\Vc e \\l [ p dx 

\JN-lJ\m\=r J JN-l<\x\<N 

< tt || Vc.||2o -^^TT" / \x\pAx<C 2 t I \x\pdx 



where C f := V2tt I|Vc L|i . . 

JV=3 



W-l<|x|<JV J N-l<\x\<N 

oo r N „ / r \ 1 / 2 

\x\—r 



Since N. / / \[P l^ c c I dcr dr < C € I / |x|pdx) , it follows that 

N=3 JN-l J\x\=r \JR 2 J 

lim / / yTjlVcel da dr = 0, and thus, there exists a sequence {Rj} /* oo such that 

?^°°J N -l J\x\=r 

[ Vpl Vc e| dcr = 0. (5.48) 



lim 



Plugging (|5~47|) and ([5^48]) into (|5T46| . we get 

(II) < 2 / y^j (j e *p* 7 e ) da; + 32 tt + 4v^(J 7 + 2tt) . 
Finally, a simple application of Holder's inequality gives 

/ \fp~ilt *P*lc)dx < || % /p||3||7e*P*7e||3/2 < / p 3/2 dx , 
JK 2 JB 2 

to conclude 

(II) < 2 y p 3/2 da; + 32 vr VS + AV5( J 7 + 2tt) . (5.49) 

Control of (III): Remind that f' s (s) < (8A)~ 1/2 := n x /2 and that 2f' s (s) -)• k a as 5 -> 0, see Proposition |2~2"1 By 
definition of c £ and G e and by symmetry of 7 



(III) = / p{x) f' s (x) x-{VG c *p) (x) dx = / / p(x) f' s (x) x 7e (x - z) ( VG * 7e * p) (z) dz dx 

JJr 2 xr 2 

(pf' 5 id*j e )(z) (VG* 7t *p) (z) dz . 

By definition of the convolution, we have 

pf 5 id*"/ e (z) = / p(z - x) f' 5 {z - x) (z - x) ~f e (x) dx 
Js. 2 

= z / p(z - x) f' s (z - x) 7e (x) da; - / p(z - x) f' s (z - x) x 7e (x) da; 

il 2 JR 2 
= Z (Pfs * 7e)(z) ~ (Pf's * id 7<0M ■ 

As a consequence, we infer 

(HI) = / [z (pf s * j e )(z) - (pf s *id le )(z)}-(VG*j e *p)(z)dz:= (111)^(111)^ (5.50) 
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By the symmetrization argument just as in (|1.13p . the first term of the right hand side of (|5.50[) reads 

(HI)! = I [ z (p(2f s - k x ) * le)(z) ■ (VG * 7e * P) (z) dz - Sttkx 
:= (III) U - 87TKA ■ 

We now control (III) 1:L using the HLS and Young inequalities, see (|4.6I) to obtain 

|(III) n |<i // \z\\(p(2f> s - Kx )* le )(z)\-^—\( lt *p)(x)\dzdx 

1 J JR 2 xR 2 \ x ~ z \ 

^ UN \(p\ 2 f's - «A| * 7e)IU/3 IIP * Te|U/3 

< ^ (IKN P|2/i " * 7.IU/3 + ||(p|2/i - * (\X\ 7e)||4/ 3 ) ||/> * 7e|U/3 

< 2C HLS (||M p|2/^ - Kx\\\l he\U/3 + WpW's - UN TelU/s) IKIU/3 , 

from which 

(III) a < Sttkx + C.W |2/i-K A |(l + |a:|)p||i. (5.51) 
To estimate the second term of the right hand side of (|5.50[) . we make again use of the HLS inequality, see 

|(III) 2 |<_ L- ff (p*| z | 7£ )( z )_^ 7£ * p )(y) d ydz 
4V2Att JJr2 xR 2 \z-y\ 

< (M7e)IU/3 Il7e * P\U/3 • 



4V2A7T 

By the Young inequality, and a direct calculation, \\p* (\z\ O'e) IU/3 < 87r || |x| II 4/3 = 8tt yfk || |x| 7IU/3- and in the 
same way ||7 e * pWa/z < 8tt ||7 e ||4/3 = 8tt — = 1 1 T 1 1 4 /3 - The positive and negative powers of e cancel, and using (|5.5ip . 
we conclude 

(III) < -8tt Ka + C t || \2f' & - nx\ (1 + \x\)p\\i + 16 ^ LS |||x| 7II4/3 H-YlU/3 

V 2A 



-Sttkx + C t || \2f' s - k x \ (1 + N)p||i + (5.52) 

v 2A 



Let us estimate this third term in a different way that will be useful later on. Using again the Young inequality, 
but this time eliminating 7 instead of p, i.e. ||7 e * PH4/3 < HpIU/3, we get 

IIP* (M7 e )ll4/3ll7 e *p|U/3 < 647r 2 Ve|||a;|7||4 / 3 ||p|| 4/ 3 . 
As a consequence, we get this other control on (III) by 

(III) < -8™ A + C e || \2f' s -Kx\(l+ \x\)ph + 16 ^ LS V~e II W 7IU/3 IHU/3 • (5.53) 

V 2A 



Control of (IV): By integrating by parts for any R > 0, we have 

xf' s (\x\ 2 )-Vp(x)dx = [ p(x)f s (\x\ 2 )x-nda- [ V • [x f s (\x\ 2 )} p(x) dx , 

J\x\=R J\x\<R 



'\x\<R J\x\=R J\x\<R 

where n denotes the outward normal to the disk Dr. Taking into account that 

~ 3 JN-lJ\x\=r 2\/2A ~ 3 JN-l<\x\<N 



E/ / \x\\f s (\x\ 2 )\p(x)d*dr<— = ]T / \x\pdx< — = / \x\pdx<^ 
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we have lim / / p{x) |/^(|a;| 2 )| |x| da dr = 0. As a consequence, there exists a sequence {Rj} /~ oo such 

w ^o° Jn-1 J\x\=r 

that lim / p{x) \f'$(\x\ )\ \x\ da = 0, and thus, we conclude 

(IV) = jj-[xf' s (\x\ 2 )]p(x)dx. (5.54) 
The desired estimates are obtained by putting together estimates (|5.49p . (|5.52l) . (|5.53[) and (|5.54[) . ■ 
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